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A  Survey  of  the  Turbulence  in  the  Marine  Surface  Layer  for  the 
Operation  of  Low-Reynolds  Number  Aircraft 

1.0  Introduction 

The  motivation  for  the  preparation  of  this  report  is  the  need  for 

low-level  wind  input,  required  for1  _the  design  and  operation  of  remotely 

piloted  aircraft  at  low-Reynolds .number^  an  at  extremely  low  altitudes 
v-e  A >  r  '■*>  n\  a.  tripod  cJoi  . 

above  the  oc£an^  Surface  winds  and  gusts  are  known  to  have  a  strong 

influence  on  the  operation  of  aircraft  near  the  surface  and  can  cause 
serious  control  problems.  In  general,  cross  winds,  wind  shear  and  gusts 
encountered  by  low-flying  aircraft  pose  significant  problems  in  their 
operation. 

However,  most  analytical  and  experimental  work  on  the  atmospheric 
boundary  layer  over  the  ocean  has  been  directed  toward  the  understanding 
of  the  processes  dealing  with  the  transport  of  mass,  momentum  and 
energy.  The  majority  of  surface  wind  measurements  has  been  collected  by 

meteorologists  who  study  these  processes,  and  has  not  been  collected 

\ 

with  the  needs  of  the  aeronautical  engineer  in  mind.  -'Although  numerous 
experimental  studies  have  been  made  in  the  surface  layer  and  boundary 
layer  over  the  ocean,  the  ability  to  apply  the  available  information  to 
the  design  and  operation  of  aircraft  is  limited. 

Notwithstanding  this  limitation,  this  report  reviews  the  current 
experimental  knowledge  of  the  marine  atmospheric  surface  layer  and 
discusses  the  nature  and  quality  of  the  reported  data  as  well  as  their 


applicability  as  input  to  the  design  and  operation  of  low-flying 
aircraft. 


2.0  General  Discussion  of  the  Problem 

For  aeronautical  applications  at  high  altitudes,  local  wind  speed 
Is  relatively  small  and  for  design  purposes  the  velocity  is  taken  as 
that  of  the  aircraft.  For  low  altitude  situations,  particularly  with 
the  speed  of  the  aircraft  being  of  the  same  order  as  the  wind  speed,  the 
relative  velocity  between  the  aircraft  and  the  wind  should  be  used.  For 
aircraft  flying  at  extremely  low  ground  speeds  of  approximately  10  m/s 
at  altitudes  of  50m  or  less,  the  atmospheric  winds  will  not  only  affect 
the  overall  operation  and  control  of  the  aircraft  but  also  affect  the 
airfoil  performance.  Under  these  conditions  the  relative  wind  velocity 
as  experienced  by  the  aircraft  will  alter  dramatically  as  the  aircraft 
moves  through  the  highly  variable  wind  field. 

The  vertical  turbulence,  cross-winds,  lateral  wind  shear  and  gusts 
encountered  by  this  type  of  aircraft  will  undoubtedly  affect  its 
performance.  However  there  does  not  seem  a  complete  understanding  how 
wind  and  gusts  enter  into  the  design,  operation  and  control  of  any 
aircraft.  Consequently,  it  is  not  totally  clear  what  information  of  the 
ground  winds  is  needed  to  aid  the  aeronautical  engineer  in  the  design  of 
these  unusual  aircraft. 

In  addition,  the  atmospheric  turbulence  can  have  a  direct  affect  in 
many  different  ways  on  the  airfoil  performance  at  chord  Reynolds  numbers 
of  the  order  of  100,000.  Under  these  conditions  angle  of  attack, 
separation  and  transition  from  the  laminar  to  the  turbulent  regime  of 
either  the  foil  boundary  layer  or  wake  have  a  pronounced  effect  on  the 
lift  and  drag  of  the  airfoil.  Especially  when  the  chord  Reynolds  number 
decreases  below  100,000  can  abrupt  performance  changes  be  expected  [11. 


It  Is  the  concensus  among  aeronautical  engineers  that  the  overall 
performance  of  the  aircraft  Is  much  affected  by  the  vertical  force 
spectrum  and  pitching  moment  spectrum  due  to  u  and  w  turbulence,  yawing 
moments  due  to  v  turbulence  and  rolling  moments  due  to  spanwlse 
variations  In  u  and  w  turbulence  (Fig.  1).  This  Is  specifically  a 
problem  for  low  Reynolds  number  flight  at  ground  velocities  of  the  same 
order  as  the  wind  speed  and  at  low  altitudes  where  there  Is  little  room 
for  error. 

Of  course,  the  fluctuations  of  the  rolling,  pitching  and  yawing 
motion  of  the  aircraft  are  primarily  caused  by  turbulent  eddies  of  a 
scale  comparable  to  the  dimension  of  the  aircraft  [3].  On  the  other 
hand  the  performance  of  the  airfoil,  controlled  by  the  boundary  layer 
flow  over  the  foil.  Is  primarily  governed  by  turbulence,  fluctuations  of 
much  smaller  scale,  of  the  order  of  the  boundary  layer  thickness. 

However  detailed  Information  of  the  turbulence  at  these  small  scales  and 
of  the  non-uniformity  of  the  gusts  over  the  dimensions  of  the  aircraft 
(cross-correlations)  Is  extremely  scarce  especially  for  the  atmospheric 
flow  over  the  ocean. 

Existing  data  obtained  over  the  ocean  are  either  from  Instruments 
mounted  on  floating  platforms  primarily  below  20m  or  from  Instrumented 
aircraft  above  20m.  Moreover,  the  marine  environment  is  not  very 
conducive  to  the  acquisition  of  multipoint  turbulence  measurements,  and 
available  turbulence  data  obtained  from  floating  platforms  are  usually 
limited  to  a  single  level.  Only  a  limited  amount  of  multipoint 
turbulence  data  Is  available  from  micrometeorologlcal  towers  located  on 
land  close  to  the  water  for  on-shore  winds. 


The  results  presented  In  this  report  are  obtained  from  many 
different  sources  dealing  with  mean  wind  and  turbulence  measurements  In 
the  atmospheric  surface  layer  over  the  ocean.  Particular  emphasis  Is 
placed  on  the  Identification  and  analysis  of  the  low-altitude  turbulence 
measurements  that  have  the  greatest  influence  on  the  operation  of  the 
low  Reynolds  number  aircraft  In  the  marine  surface  layer. 

3.0  General  Discussion  of  the  Atmospheric  Boundary  Layer 

The  atmospheric  or  planetary  boundary  layer  (A.B.L.  or  P.B.L.)  is 
the  result  of  the  Interaction  of  the  atmospheric  flow  over  the 
underlying  land  or  sea  surface.  This  layer  Is  characterized  by  a 
turbulent  transfer  of  momentum,  heat  and  mass  (water  vapor)  and  their 
associated  gradients.  Non-uniform  boundary  conditions, 
nonstatlonarltles,  variable  stability  and  thermal  conditions  as  well  as 
Coriolis  effects  are  conditions  which  are  typical  for  the  A.B.L.  but  not 
common  to  wind  tunnel  or  laboratory  boundary  layers.  Moreover,  under 
extreme  Inversion  and  convective  conditions,  the  latter  Including  the 
passage  of  fronts,  squall  lines,  thunderstorms  and  other  mesoscale 
phenomena,  the  surface  flow  cannot  be  classified  as  typical  boundary 
layer  flow. 

During  daytime  under  generally  convective  conditions  the  boundary 
layer  depth  is  assumed  to  correspond  to  the  altitude  of  the  lowest 
inversion,  z^,  typically  about  1000  m  [41.  Away  from  the  surface  the 
large  scale  buoyancy-generated  turbulence  promotes  the  transport  and 
mixing  of  heat  and  momentum.  As  a  result  of  this  vigorous  mixing  the 
mean  profiles  In  the  upper  90%  of  the  A.B.L.  below  the  inversion  are 
almost  uniform  (mixed  layer),  while  above  the  Inversion  base  the 


velocity  and  potential  temperature  Increase  with  height.  Another  factor 
which  complicates  the  nature  of  the  A.B.L.  Is  the  downward  flux  of 
momentum  and  heat  due  to  entrainment  Into  the  convective  boundary  layer. 

Ourlng  nighttime,  thermal  conditions  are  often  stable  near  the 
surface  and  the  depth  of  the  surface  inversion  Is  considered  to 
correspond  to  the  A.B.L.  (typically  about  100  m  [4]).  Under  these 
conditions  the  surface  cools  down  resulting  in  a  stable  temperature 
stratification  just  above  the  surface  where  the  buoyancy  forces  will 
Inhibit  the  vertical  motion  and  suppress  the  turbulent  mixing.  Above 
the  A.B.L.  the  flow  Is  considered  frictionless  and  primarily  governed  by 
pressure  forces,  Coriolis  forces  and  horizontal  temperature  gradients. 

The  surface  layer,  approximately  the  lowest  10%  of  the  A.B.L.,  Is 
characterized  by  the  mechanically  produced  turbulence  from  the  surface 
roughness  and  friction,  with  nearly  constant  vertical  fluxes  of 
momentum,  heat  and  mass.  The  height  of  this  layer  may  vary  over  a  wide 
range  but  Is  typically  of  the  order  of  50  m  [4].  In  the  lower  most  part 
of  the  surface  layer,  commonly  referred  to  as  the  roughness  layer,  the 
flow  Is  affected  by  the  Individual  roughness  elements  and  is 
nonhomogeneous  and  three-dimensional  in  nature.  Several  roughness 
heights  above  the  mean  surface,  this  three-dimensionality  of  the  flow 
disappears  and  the  flow  Is  usually  considered  to  be  horizontally 
homogeneous.  In  the  marine  roughness  layer  various  properties  of  the 
sea  surface,  primarily  the  travelling  surface  waves  affect  the  flow. 
However,  the  observation  of  these  wave  effects  have  been  hindered  by  the 
motion  of  the  instrument  platform  induced  by  the  wind  and  waves  as  well 
as  by  the  presence  of  the  platforms  (towers,  buoys,  ships)  in  the  flow 


There  exists  a  strong  Interaction  between  the  oceanic  and 
atmospheric  motions.  The  atmospheric  flow  Is  responsible  for  the  ocean 
currents,  bringing  layers  of  water  of  different  temperature  to  the 
surface,  which  In  turn  affects  the  air  temperature  and  the  water  vapor 
content  and  thus  affects  the  atmospheric  motion.  The  trade  winds  drive 
the  ocean  currents  north  and  south  of  the  equator,  while  a  combination 
of  high-pressure  and  low-pressure  cells  cause  the  warm  sea  currents  like 
the  Gulf  Stream  and  cold  sea  currents. 

4.0  The  Governing  Equations 

The  basic  conservation  laws  of  mass  momentum  and  energy  are  used 
for  the  description  of  the  motion  of  the  atmosphere.  The  momentum 
equation  should  Include  body  and  buoyancy  forces,  Coriolis  forces  In 

addition  to  the  pressure  and  viscous  forces.  The  use  of  this  set  of 

equations  may  be  greatly  simplified  by  adopting  the  Bousslnesq 
approximation.  In  this  approximation  the  atmosphere  is  assumed  to  be  a 
perfect  gas,  while  pressure,  density  and  temperature  vary  only  slightly 
from  the  neutrally  stable  atmosphere  (hydrostatic  equilibrium). 
Temperature  variations  are  considered  to  be  sufficiently  small  so  that 
viscosity  and  heat  conductivity  can  be  considered  constant.  Also  It  Is 
assumed  that  the  velocity  of  the  flow  is  sufficiently  small  so  that 
variations  in  density  due  to  variations  in  dynamic  pressure  are  small. 
However,  the  variations  In  density  due  to  variations  in  temperature,  as 
the  fluid  parcels  change  elevation  and  adapt  themselves  to  the  new 
environment,  cannot  be  Ignored.  These  variations  in  temperature  and 
density  produce  a  non-equilibrium  between  buoyancy  and  gravitational 
forces.  The  net  amount  of  work  done  by  these  forces  may  be  positive 
(turbulent  production)  or  negative  (turbulent  suppression).  Although 


generally  speaking  the  A.B.L.  flow  Is  purely  turbulent.  In  stable  air 
non-turbulent  gravity  waves  will  also  produce  velocity  fluctuations. 


Considering  the  flow  to  be  turbulent  and  with  Reynolds  averaging  of 
the  Instantaneous  equations,  the  continuity  equation  for  the  mean  flow 
becomes 

3U, 

— 1  =  0  (Incompressible  flow)  (1) 

9X  ^ 


Here  capital  letters  denote  averages  and  lower  case  letters  denote 
fluctuations.  Similarly,  the  mean  momentum  equation  becomes 


1  aP'  ,  1 
°o  3X1  0 
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Here  the  subscript  o  indicates  properties  associated  with  the  reference 
state,  and  the  primes  Indicate  deviations  from  this  reference  state. 

This  form  of  the  momentum  equation  Includes  the  Reynolds  stress  term  and 
the  resultant  between  the  gravitational  and  buoyancy  force.  The  thermal 
energy  equation  becomes 


0 T 1  _  I  3  r  v  3 T 1  _  m  /  \ 
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which  Include  the  conduction  (usually  negligible),  the  turbulent  heat 
flux  and  the  heat  flux  due  to  radiation.  Because  of  the  presence  of  the 
turbulent  flux  terms  this  set  of  equations  can  no  longer  be  considered 
to  be  a  closed  set. 


So  far  the  effect  of  moisture  fluctuations  on  the  density  has  been 
neglected,  which  might  be  quite  reasonable  for  A.B.L.  flow  over  land. 
However  over  water  and  specifically  over  warm  water  (tropical  oceans) 
the  effect  of  moisture  is  not  negligible.  Under  these  conditions  the 
variations  in  specific  humidity  can  be  Incorporated  in  the  governing 


equations  by  replacing  the  temperature  and  heat  flux  by  virtual 
temperature  and  virtual  heat  flux  (51.  The  virtual  temperature  is 
defined  as  the  temperature  dry  air  would  have  In  order  to  obtain  the 
observed  density  and  pressure  of  the  moist  air  In  which  the  observations 
are  being  made. 

Equations  for  the  turbulent  terms  and  specifically  the  turbulent 
kinetic  energy  e  *  ^  u^  can  also  be  obtained.  With  Reynolds  averaging 
and  the  neglectlon  of  the  small  viscous  term  and  the  assumption  of 
horizontal  homogeneity,  the  equation  of  mean  turbulent  energy,  E  becomes 

1  2  3  4 

Each  of  the  terms  on  the  rlghthand  side  has  the  following  physical 
Interpretation: 

1.  The  mechanical  production  of  turbulence. 

2.  The  rate  of  turbulence  production  ew  >  0  for  unstable 
atmosphere,  or  the  rate  of  suppression  of  turbulent  energy 
ew  <  0  for  a  stably  stratified  atmosphere. 

3.  The  viscous  dissipation. 

4.  The  diffusion  of  turbulent  energy. 

The  ratio  of  terms  2  and  1  represents  a  measure  of  the  buoyancy 
production  or  suppression  of  turbulence  relative  to  the  mechanically 
produced  turbulence,  and  is  referred  to  the  flux  Richardson  number 


T  —  dU 
0  uw^ 


(5) 


For  vanishing  R^  the  atmospheric  flow  Is  considered  to  be  neutrally 
stable,  for  Rf  <  0  with  ew  >  0  (upward  heat  flux)  the  atmosphere  is 


thermally  unstable  and  for  Rf  >  0  with  ew  <  0  (downward  heat  flux)  the 
atmosphere  Is  thermally  stable. 

With  the  Introduction  of  turbulent  dlffuslvltles  or  transport 
coefficients  and  assumption  that  these  coefficients  are  equal,  the 
gradient  Richardson  number  can  be  defined  as 

R1  =  f  (6) 

'o  (dU/dzr 

Since  no  detailed  solution  of  the  turbulent  flow  can  be  expected  from 
the  above  set  of  equations,  they  can  only  provide  a  framework  for  the 
Interpretation  and  description  of  observations  and  a  basis  for  the 
modelling  of  the  flow. 

5.0  Similarity  Theories 

The  description  of  turbulent  boundary  layers  has  been  based 
primarily  similarity  arguments.  Using  this  approach,  the  statistical 
flow  parameters  are  described  In  terms  of  empirical  functions  which 
require  a  number  of  basic  parameters  which  control  the  flow  In  different 
parts  of  the  boundary  layer.  Specific  empirical  relations  have  been 
derived  from  experimental  data,  by  organizing  the  data  in  accordance 
with  the  established  turbulent  boundary- layer  model. 

5.1  Wind  Tunnel  Boundary  Layers 

5.1.1  Smooth  Boundary 

Based  on  dimensional  analysis  and  experimental  results  from 
turbulent  boundary  layers  In  the  laboratory  two  basic  laws  emerged.  The 
"law  of  the  wall"  describing  the  mean  flow  near: 

a.  smooth  wall  U/U*  *  f(-^-jp) 


(7) 


b.  rough  wall  U/UQ  *  f(z/k). 


and  the  "velocity  defect  law"  governing  the  mean  flow  in  the  outer  layer 
for  either  rough  or  smooth  walls 

(U  -  UQ)/U*  -  g(z/6).  (9) 

Millikan  showed  that  the  existence  of  an  overlap  region  of  the  two  laws 
requires  a  logarithmic  velocity  profile  without  introducing  eddy 
viscosity,  mixing  length  and  without  assuming  the  existance  of  a 
constant  stress  layer. 

Using  the  method  of  asymptotic  matching  Tennekes  and  Lumley  [61 

showed  the  existence  of  the  logarithmic  velocity  profile  for  smooth 

flat-plate  boundary  layers  provided  the  Reynolds  number  is  large 
* 

(Ugfi/v  -  <»).  Under  these  conditions  the  velocity  profiles  are  self¬ 
preserving  and  the  normalized  velocity  defect,  (U  -  UQ)/U*,  and  the 
normalized  turbulent  stresses  are  universal  functions  of  z/d, 

* 

Independent  of  downstream  distance  or  any  other  parameter.  As  UQd/v  or 

s/(v/U*)  ♦  ®,  the  law  of  the  wall  can  be  matched  with  the  velocity 

★ 

defect  law  in  the  overlap  region  where  simultaneously  z  »  v/UQ  and  z 
«  6.  In  this  region  the  nondimensional  velocity  gradients  obtained 
from  either  the  law  of  the  wall  or  the  velocity  defect  law  must  be  equal 
and  also  independent  of  any  other  variables  and  therefore  equal  to  a 
universal  constant.  Integration  of  both  velocity  gradients  and 
evaluation  of  the  constants  from  experimental  results  lead  to 

U/U*  =  2.5  in  (U*z/v)  +  4.9  (10) 


(U  -  UJ/U„  =  2.5  tn  (z/6)  -  2.5 


[i] 


5.1.2  The  Effect  of  Roughness 

The  effect  of  wall  roughness  on  the  flow  In  the  boundary  layer  can 

only  be  treated  In  general  terms  because  of  the  variation  In  geometrical 

shapes  and  distribution  of  the  roughness  elements.  The  roughness  does 

not  directly  affect  the  velocity  defect  law  except  through  the  value  of 

the  wall  shear  or  friction  velocity,  U*,  provided  the  size  of  the 

uniform  roughness  elements,  k,  is  much  smaller  than  the  boundary  layer 

★ 

thickness,  s.  The  friction  velocity,  UQ,  can  no  longer  be  considered  to 
be  a  local  value  but  instead  must  be  interpreted  as  an  average  value 
over  an  area  which  includes  a  sufficient  number  of  roughness  elements. 

On  the  other  hand,  the  velocity  distribution  close  to  the  boundary  Is 
affected  directly  by  the  roughness,  so  that  the  law  of  the  wall  may  be 
written  as 

U/Uo  =  f<f  ■  02) 

In  the  overlap  region  where  z  »  k  and  z  «  5,  the  nondimensional 
velocity  gradients  must  be  equal  to  the  same  universal  constant  as  in 
the  smooth-wall  case.  Integration  leads  to 

*  U*k 

U/UQ  =  2.5  an  (z/k)  +  F(-S-)  (13) 

where  F  is  not  a  universal  constant  but  for  a  given  roughness  size  also 

depends  on  the  shape,  distribution  and  density  of  the  roughness  elements 
★ 

even  If  k  »  v/UQ,  fully  rough  wall.  Also  note  that  the  logarithmic 
profile  is  valid  only  in  the  overlap  region  where  z  »  k.  In  the  region 
where  z  and  k  are  of  the  same  order  (roughness  layer)  the  above 
criterion  is  not  satisfied,  and  the  logarithmic  profile  should  not  be 
expected  to  apply. 


5.2  Atmospheric  Boundary  Layer 

5.2.1  Neutral  Atmospheric  Boundary  Layer 


Consider  the  flow  in  the  A.B.L.  as  being  steady  and  horizontally 

homogeneous  without  the  effect  of  any  buoyancy  forces  and  over  a  flat 

homogeneous  surface  whose  roughness  can  be  represented  by  a  roughness 

length-scale,  zQ.  Above  the  A.B.L. ,  the  pressure  gradient  force  must  be 

balanced  with  the  horizontal  component  of  the  Coriolis  force  and 

fC  x  G  =  -  i  vp  (14) 

where  f  is  the  Coriolis  parameter  f  =  2<dSlnx,  with  u  being  the  earth's 

angular  velocity  and  x  the  geographic  latitude  and  G  being  the 

geostrophic  wind.  The  vertical  component  of  the  Coriolis  force  is 

negligible  In  comparison  with  gravity.  The  external  parameters  which 

control  the  flow  are  f,  G  and  zQ.  Other  flow  parameters  which  enter  are 

the  friction  velocity  U*  and  the  height  of  the  neutral  A.B.L.,  h.  In 

the  absence  of  buoyancy  forces,  the  turbulence  is  purely  of  mechanical 

* 

origin  and  the  friction  velocity,  UQ,  Is  the  only  characteristic 
velocity  scale  for  the  turbulence  (7].  The  height  h,  which  is 

k 

proportionally  to  UQ/f,  and  zQ  must  be  the  characteristic  length  scales 

for  the  flow  in  the  outer  and  inner  layer  respectively.  In  the  limit  as 
★ 

U0/fzQ  ♦  ®  the  two  similarity  laws  one  for  the  surface  layer  and  one  for 
the  outer  layer  must  match  in  the  overlap  region  where  both  laws  are 
valid  simultaneously.  In  this  region  with  z/zQ  ♦  »  and  zf/U*  ♦  0  the 
nondimensionalized  velocity  gradients  should  be  equal  to  the  same 
universal  constant  as  before.  Integration  leads  again  to  the  well  known 
logarithmic  velocity  profile 

VJm  (z/zj  +  C  (15) 


Although  at  the  present  time  some  controversy  exists  over  the  numerical 

value  of  the  value  of  K  (the  von-Karman  constant),  a  value  of  K  =  0.4 

based  on  wind  tunnel  experiments  Is  most  commonly  accepted.  As  before  C 

Is  not  a  universal  constant  but  must  depend  on  the  nature  of  the 

★ 

roughness  and  the  roughness  Reynolds  number  Uqzq/v.  The  constant  of 
Integration,  C,  is  often  Incorporated  into  zQ  in  the  above  expression 

U/U*  =  i  an  (z/zQ)  (16) 

For  any  natural  roughness  it  Is  impossible  to  assign  a  single  numerical 
value  for  z0,  without  taking  Into  account  the  variation  in  roughness 
size,  shape,  distribution  and  density.  For  all  practical  situations  a 
"profile  roughness  length"  can  be  obtained  from  velocity  observations 
only  In  the  range  for  which  (16)  Is  expected  to  be  valid,  or  for  z  » 
zQ.  Also  the  friction  velocity  U*  is  certainly  not  a  local  value  but  is 
based  on  the  wall  shear  stress  rQ,  averaged  over  a  sufficiently  large 
upstream  fetch.  Using  this  derivation  of  the  logarithmic  wind  profile 
in  the  neutral  A.B.L.  does  not  require  the  assumption  of  constant  stress 
layer. 

5.2.2  Non-Neutral  Atmospheric  Boundary  Layer 

So  far  the  discussion  has  been  limited  to  a  neutral  atmospheric 
boundary  layer  and  laboratory  boundary  layers  with  adiabatic  temperature 
profiles  and  therefore  without  vertical  heat  flux.  However  the  A.B.L. 
nearly  always  has  some  vertical  heat  flux  either  in  the  form  of  sensible 
or  latent  heat.  In  addition  to  the  mechanically  produced  turbulence, 
turbulent  energy  Is  either  produced  or  suppressed  depending  on  whether 
the  heat  flux  is  upward  (unstable)  or  downward  (stable).  Consequently 


the  mean  flow  as  well  as  the  turbulence  in  the  diabatic  surface  layer 
are  different  from  those  in  the  neutral  surface  layer. 

Stable  conditions  are  usually  encountered  during  night  or  when 
warmer  air  flows  over  a  colder  water  surface  and  are  characterized  by  a 
suppression  of  the  turbulence.  Under  extreme  stable  conditions,  the  air 
in  the  different  layers  becomes  uncoupled  as  a  result  of  reduced 
mechanical  mixing,  and  jet-like  flows  either  near  the  surface  or  at 
higher  elevations  may  be  observed.  Additional  complications  may  occur 
when  turbulence  and  gravity  waves  co-exist.  Under  these  conditions  the 
flow  seldom  reaches  a  state  of  quasi-steadiness,  because  conditions  vary 
rather  rapidly,  and  consequently  no  simply  boundary- layer  model  is 
available  for  the  description  of  the  flow. 

With  daytime  heating  of  the  earth's  surface  resulting  in  an  upward 

heat  flux,  the  mechanical  turbulence  production  is  augmented  with 

buoyant  energy  production.  Near  the  surface  the  mechanical  energy 

production  dominates  but  decreases  rapidly  with  height,  while  on  the 

other  hand  the  buoyant  production  is  almost  constant  through  the  surface 

layer.  Outside  the  surface  layer,  the  unstable  or  convective  A.B.L. 

undergoes  strong  vertical  mixing  due  to  upward  buoyancy  and  the 

production  of  large-scale  convective  turbulence.  In  addition,  downward 

fluxes  of  momentum  and  heat  due  to  entrainment  into  the  convective 

boundary  layer  across  the  overlying  inversion  complicate  matters  a  great 

deal.  Near  the  surface  the  turbulent  fluxes  exist  as  a  result  of 

gradients,  however  in  the  mixed  layer  the  fluxes  are  maintained  by 

buoyancy  effects  and  entrainment.  No  single  internal  parameter  as  the 

★ 

friction  velocity,  UQ,  is  available  to  scale  the  mean  and  turbulent  flow 
in  both  layers  and  consequently  the  method  of  asymptotic  matching  cannot 


be  used.  Instead  we  have  to  resort  to  separate  similarity  theories  in 
the  surface  and  the  outer  layer. 

5.2.3  Surface  Layer  Similarity 

Monin  and  Obukhov  (8]  proposed  that  the  flow  (mean  or  turbulence) 

in  the  surface  layer  is  uniquely  determined  by  the  surface 

★ 

stress,  UQ,  the  kinematic  surface  heat  flux  Qq  =  Ho/pCp,  the  buoyancy 
parameter  g/TQ  and  the  height  z.  With  this  approach,  flow  quantities 
are  not  directly  dependent  on  the  roughness  but  are  influenced  by  it 
indirectly  through  the  magnitude  of  the  momentum  flux  or  U*.  Then  by 
dimensional  analysis,  properly  nond i mens i ona 1 i zed  flow  variables  should 
be  universal  functions  of  the  stability  parameter  z/L  where  L  is  the 
buoyancy  length  scale  known  as  the  Monin-Obukhov  stability  length. 


Accordingly,  the  dimensionless  wind  shear  should  be  a  universal  function 
of  z/L,  which  can  be  integrated  to  give 

U/U*  =  j?  Un  z/z0  '  *(Z/L)1  (18) 

where  zQ  is  introduced  as  a  constant  of  integration  and  can  only  be 
obtained  by  extrapolation  of  velocity  observations  above  the  roughness 
layer  down  to  the  axis  where  U  =  0.  Here,  the  function  i>  is  an 
empirically  derived  universal  function  of  z/L  with  the  limit  -  0  for 
the  aproach  of  neutral  stability  for  which  L  -  ®  and  z/L  -  0. 


5.2.4  Outer  Layer  Similarity 

In  the  outer  layer  the  mean  and  turbulent  flow,  besides  depending 
on  the  height  z  and  on  the  surface  fluxes  of  heat  and  momentum,  also 
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depend  on  the  boundary  layer  height,  h,  the  Coriolis  parameter,  f, 
entrainment  and  wind  shear  at  the  top  of  the  boundary  layer.  The  form 
of  the  universal  functions  describing  the  dependence  of  the  flow 
characteristics  on  these  parameters  can  only  be  obtained  from 
observations  and  numerical  modelling.  However,  most  of  the  mean  wind 
and  turbulence  observations  have  been  limited  primarily  to  the  lower  100 
m  over  land  and  to  the  lower  20  m  over  water.  Simultaneous  observations 
of  the  mean  flow  and  turbulence  throughout  the  depth  of  the  A.B.L.  are 
extremely  sparse  and  limited  to  observations  from  instrumented  airplanes 
and  from  tethered  balloons  either  with  or  without  surface  layer 
observations.  Moreover,  it  is  practically  impossible  to  study  the 
effect  of  only  one  of  the  above  parameters  on  the  flow  characteristics 
while  holding  the  remaining  parameters  constant.  Most  of  the  existing 
similarity  schemes  are  restricted  by  Including  only  one  or  two  of  the 
above  parameters  in  the  analysis. 

For  a  far  more  comprehensive  treatment  of  the  atmospheric  boundary 
layer  and  atmospheric  turbulence  the  reader  is  advised  to  consult 
References  4,  5,  9,  10  and  11. 

6.0  Discussion  of  Experimental  Procedures 
6.1  General  Observations 

Recent  Improvements  in  the  quality  of  data  acquisition,  data 
handling  and  data  processing  systems  have  resulted  in  an  increased 
number  of  experimental  research  programs  dealing  with  the  structure  and 
dynamics  of  the  turbulence  in  the  A.B.L.  Most  of  these  investigations 
have  been  performed  over  land,  using  meteorological  towers  instrumented 
at  several  levels  for  observations  generally  at  elevations  below  200m, 
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and  Instrumented  airplanes  as  well  as  instrumentation  suspended  from 
balloon  tethers  for  observations  at  higher  elevations.  These  kind  of 
observations  over  water  and  specifically  over  the  open  ocean,  where 
conditions  are  not  amenable  to  experimental  Investigations,  are  much 
fewer  In  numbers. 

Difficulties  In  obtaining  quality  measurements  of  the  turbulence 
over  water  are  generally  compounded  by  the  contradictory  requirements 
that  the  anemometers  should  be  sensitive  enough  to  measure  the  high- 
frequency  fluctuations  and  at  the  same  time  be  rugged  enough  to 
withstand  the  extreme  marine  environment.  Because  of  this  problem  wind 
observations  over  water  can  generally  be  classified  into  two  groups: 
those  based  on  mean  flow  observations  from  which  certain  turbulence 
characteristics  associated  with  the  mixing  processes  are  being  deduced, 
and  those  which  measure  mean  wind  and  turbulence  simultaneously.  Only 
data  sets  belonging  to  the  latter  category  have  been  used  for  analysis 
In  this  report.  In  addition,  observations  over  water  are  often  limited 
to  one  level,  frequently  at  a  height  of  about  10m  or  lower.  Wind 
measurements  acquired  so  close  to  the  water  surface  may  have  been  taken 
in  the  roughness  layer,  especially  at  the  higher  wind  velocities  for 
which  wave  amplitudes  may  be  lm  or  more.  Consequently,  data  acquired  at 
these  levels  vary  considerably  and  are  not  only  Influenced  by  wave 
height  but  also  by  the  direction  and  the  speed  at  which  the  waves  are 
being  propagated  [12,13,14  and  15).  It  has  been  observed  that  the  wave 
Influence  on  mean  wind  and  turbulence  occurs  to  a  height  equivalent  to 
at  least  three  wave  heights  [14j. 

Only  with  many  hours  of  observations  and  averaging  of  the  results 
can  a  representative  "average"  picture  of  the  flow  over  water  be 


obtained.  Much  of  the  observed  scatter  of  the  data  may  not  be  random 
but  may  be  due  to  variation  in  thermal  stability  or  due  to  varying  wind 
and  wave  conditions.  The  water  surface  roughness  and  hence  the 
turbulence  characteristics  above  it  seem  to  depend  on  surface-wave 
parameters  which  by  themselves  are  the  result  of  local  wind  conditions 
both  past  or  present,  but  are  also  dependent  on  the  location  of  the 
observations  with  respect  to  the  large  weather  patterns. 

Many  of  the  mean  wind  and  specifically  turbulence  observations 
suffer  from  experimental  error  such  as  contamination  by  water  spray, 
inadequate  calibration  control,  platform  oscillation,  platform  motion 
and  flow  distortion  due  to  platform  blockage.  Results  from  a  carefully 
executed  comparison  experiment  between  two  different  sonic  anemometer 
systems  show  systematic  differences  in  measured  turbulence  quantities  by 
as  much  as  40%  (16).  Another  comparison  experiment  with  three  different 
type  of  sensors  (gust  anemometer,  acoustic  anemometer  and  a  propeller 
anemometer  system)  show  variations  of  at  least  10%  even  after 
corrections  for  response  limitation  and  for  calibration  differences  were 
applied  [17].  Despite  continued  improvements  in  sensor  design  and  in 
data  acquisition  and  data  handling  systems,  one  must  expect  remote 
observations  of  the  turbulence  in  the  atmospheric  surface  layer  over  the 
ocean  to  exhibit  considerable  scatter. 

A  good  deal  of  the  scatter  may  be  of  a  strictly  statistical  nature, 
meaning  that  the  turbulence  characteristics  in  the  A.B.L.  are  by  no 
means  constant.  Time  plots  of  the  turbulent  velocity  component  show 
relatively  large  fluctuations  in  relation  to  the  average  value. 
Consequently  with  finite-time  averaging,  the  presence  of  low-frequency 
fluctuations  will  lead  to  a  considerable  scatter  of  the  averaged 


quantities.  In  addition,  turbulence  quantities  measured  at  a  given 
location  can  not  be  Interpreted  as  being  an  accurate  estimate  of  what 
the  turbulent  flow  Is  at  any  other  location  in  a  large  area  of  several 
hundred  kilometers  in  dimension.  Therefore,  the  turbulence  measurements 
at  a  given  location  at  a  given  Instant  can  only  be  expected  to  be  a 
small  sample  out  of  a  set  of  infinite  samples  which  makes  up  the  process 
of  atmospheric  turbulence. 

Actually,  the  conditions  of  steady  state  and  spatial  homogeneity, 
which  are  necessary  requirements  for  similarity  theory,  are  seldomly 
satisfied  for  100  percent.  Consequently,  the  data  presented  in  this 
report  can  only  be  considered  as  estimates  based  on  a  limited  number  of 
observations.  In  view  of  this  and  with  all  the  other  uncertainties  of 
Instrumentation  performance,  data  analysis,  platform  motion,  etc.,  it 
would  be  Impossible  to  expect  the  results  presented  in  this  report  to 
show  perfect  uniformity. 

For  an  excellent  review  of  the  different  methods  and  their 
limitations  for  the  measurement  of  turbulence  fluxes  In  the  atmospheric 
surface  layer  over  water,  the  reader  should  consult  Reference  18. 

6.2  Turbulence  Stress  and  Friction  Velocity 

An  abundance  of  turbulent  flux  coefficients  are  presented  in  the 
open  literature,  however,  they  are  usually  referenced  to  a  height  of 
10m.  These  results  are  primarily  acquired  from  instrumentation  mounted 
on  masts  on  floating  platforms,  on  fixed  towers  in  relatively  shallow 
coastal  waters,  lakes  or  on  beaches,  or  from  instrumentation  mounted  on 
ships  or  airplanes  over  the  open  ocean.  For  most  of  these 
investigations  the  turbulence  observations  are  limited  to  the  turbulent 
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stress  and  heat  fluxes  (sensible  and  or  latent).  There  are  several  ways 
to  measure  the  turbulent  fluxes  of  momentum,  heat  and  moisture  In  the 
marine  surface  layer.  Concentrating  on  the  momentum  flux  or  turbulent 
stress  (-ouw),  or  the  friction  velocity,  U*  =  ^/-uw,  one  can  obtain 
values  of  these  parameters  as  follows. 

The  first  method  Is  the  profile  or  gradient  technique  in  which  the 
profile  friction  velocity,  U*,  Is  estimated  from  mean  velocity 
measurements  at  several  levels  above  the  roughness  layer,  but  still  in 
the  overlap  region  where  the  logarithmic  velocity  law  (18)  with  the 
appropriate  stability  correction  is  valid.  The  recommended  stability 
corrections  as  derived  from  land  observations  are  i|>  =  -5  z/L  for  stable 
air  [11],  and  are  given  by  Paulson  [19]  for  unstable  conditions. 

However  the  stability  parameter  z/L  itself  requires  the  knowledge  of  the 
friction  velocity  U*.  8y  using  additional  empiricism,  z/L  can  be 
estimated  from  the  gradient  Richardson  number  [111,  the  latter  requiring 
the  measurement  of  the  mean  velocity  and  the  virtual  potential 
temperature  at  two  levels  above  the  roughness  layer.  In  addition  to  the 
stability  correction  to  the  profile  there  is  also  an  uncertainty  about 
the  value  of  the  vonKarman's  constant  (0.35  <  K  <  0.41),  which 
obviously  will  affect  the  values  of  the  profile  friction  velocity. 

Blanc  [20]  provides  estimates  for  the  uncertainty  of  the  profile-flux 
measurements  by  scaling  of  the  results  from  one  experiment  to  another. 
This  analysis  showed  that  the  variation  for  typical  values  of  the 
turbulent  stress  and  the  roughness  length  fell  in  the  range  from  10  to 
40%  and  25  to  100%  respectively.  To  say  the  least,  the  profile  method 
is  quite  empirical  and  both  Krugermeyer  et  al.  [14]  and  Wieringa  [21] 
conclude  that  the  aerodynamic  stress  can  not  be  estimated  with 
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sufficient  accuracy  from  wind  profile  measurements  close  to  the  water 
surface. 


The  second  method  for  obtaining  the  friction  velocity,  U*,  is  by 
the  dissipation  technique,  for  which  the  measurement  of  the  downwind 
velocity  spectrum  at  a  single  elevation  is  required.  Successful 
application  of  this  method  requires  delicate  instrumentation  with  very 
small  sensors  and  a  high  frequency  response  up  to  2000  Hz.  The  method 
Is  based  on  the  assumption  that  the  production  of  mechanical  energy  Is 
equal  to  the  dissipation  for  near  neutral  stability  conditions  or 
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With  the  assumption  of  the  existence  of  the  logarithmic  velocity 
profile.  In  differentail  form 

dU  uo 

Sz  ■  Ki  <20> 

and  assuming  that  the  turbulent  stress  is  constant  with  height  and  equal 

to  the  profile  stress 

-uw  =*  (U*)2  =  (U*)2  (21) 


one  obtains 


(U*)z  *  (Kez)2/3  (22) 

The  value  of  the  dissipation,  e,  can  be  obtained  from  the  measured 
velocity  spectrum  which  in  the  inertial  subrange  should  follow 

*u(k)  =  <  e2/V5/3,  (23) 

where  *  is  the  Kolmogoroff  constant  (<  =  0.5)  and  k  =  ,  is  the 

radian  wave  number.  The  values  of  both  the  vonKarman  constant,  K,  and 
the  Kolmogoroff  constant,  *,  are  still  somewhat  controversial.  In  the 
equilibrium  subrange  the  turbulence  is  assumed  to  be  locally  isotropic 
which  requires  the  value  of  the  ratio  *w(k)/$u(k)  to  be  4/3.  However 
measurements  8m  above  the  water  surface,  by  Schmitt  et  al.  (221  show 


values  of  this  ratio  Instead  to  be  close  to  unity.  Indicating  a 
condition  of  anisotropy  of  the  turbulence  in  the  inertial  subrange. 

Also  under  non-neutral  stability  conditions  expression  (22)  must  be 
modified  requiring  an  estimate  of  the  stability  parameter  z/L  as 
discussed  before.  Besides  requiring  delicate  instrumentation  for  the 
measurement  of  the  velocity  spectra  at  sufficiently  high  frequencies, 
this  method  Is  likely  to  suffer  from  the  fact  that  basic  assumptions  may 
not  be  satisfied. 

The  third  and  possibly  the  most  frequently  used  method  for 
measuring  the  momentum  flux  is  the  eddy-correlation  technique  which 
requires  the  direct  measurement  of  the  covariance  of  the  vertical  and 
downwind  velocity  components.  This  method  Is  also  used  for  the 
measurement  of  the  variances  of  all  three  velocity  components,  requiring 
delicate  sensors  capable  of  measuring  these  components  up  to  a  frequency 
of  at  least  10  Hz.  The  variances  and  covariances  can  be  obtained  by 
multiplication  of  the  instantaneous  values  and  consequent  time 
Integration  or  may  be  obtained  by  integration  of  the  spectra  and 
cospectra  respectively.  Both  methods  suffer  from  the  same  problem, 
namely  that  the  contribution  of  the  low-frequency  fluctuations  to  the 
variances  and  covariances  Is  not  well  established.  For  the  majority  of 
the  reported  observations,  only  fluctuations  between  0.001  Hz  and  1  Hz 
are  Included  In  the  evaluation  of  these  quantities.  Inspite  of  the 
delicate  sensor  requirement,  limited  frequency  band-width,  large  amount 
of  raw  data  acquisition  as  well  as  the  sensitivity  of  the  momentum  flux 
to  sensor  alignment,  most  research  programs  dealing  with  turbulence 
measurement  have  used  this  method.  Certainly  the  eddy  correlation 
method  for  measuring  the  momentum  flux  requires  fewer  assumptions  than 


22 


the  profile  and  dissipation  methods.  Consequently,  the  majority  of  the 
turbulence  data  Included  In  this  report  used  the  eddy-correlation  method 
for  analysis,  and  were  obtained  primarily  by  either  sonic  anemometer  or 
thrust  anemometer,  although  In  some  of  the  referenced  research  programs, 
propeller  or  cup  anemometers  were  used.  In  the  case  fluxes  obtained 
from  two  different  methods  were  available,  the  preference  was  given  to 
those  measured  by  the  eddy-correlation  technique. 

7.0  Discussion  of  Experimental  Results 

The  flow  In  the  atmospheric  boundary  layer  is  usually  considered  to 
be  horizontally  homogeneous  and  steady.  Ignoring  the  effects  of 
entrainment,  cloud  cover,  Coriolis  acceleration  and  barocllnlty,  the  ABL 
can  be  divided  into  2  distinct  regions  each  with  Its  own  governing 
parameters.  In  the  surface  layer  the  mean  flow  and  turbulence  are 
governed  by  the  surface  stress,  tq,  the  elevation  above  the  surface,  the 
buoyancy  parameter  g/Ty,  surface  heat  flux  HQ  and  the  flow  properties 
pq  and  cp.  Consequently,  flow  parameters  nondlmensionalized  with  the 
shear  velocity,  U*  =  /xo/ pQ,  should  be  universal  functions  of  the 
stability  parameter,  z/L,  where  the  Monin-Obukhov  length,  L 
-  -U*3/K(g/Tv)(H0/o0cp). 

In  the  outer  layer  the  flow  is  often  no  longer  influenced  by  the 
shear  but  is  usually  driven  by  buoyancy.  Under  strong  convective 
conditions  when  the  buoyancy  forces  promote  strong  mixing,  a  so-called 
mixed  layer  Is  formed  where  the  turbulence  is  governed  by  z^,  the  height 
of  the  mixed  layer,  g/Ty,  HQ,  pQ  and  Cp.  The  velocity  scale  for  this 
regime  is  the  convective  velocity,  W*  =  l(g/Tv)  ;v°oV  zi|1/3- 
Consequently,  turbulence  parameters  nondlmensionalized  with  W*  are 


expected  to  be  universal  functions  of  z/z^  only.  Over  land  where  the 
heat  fluxes  are  usually  much  larger  than  over  the  oceans,  the  surface 
layer  and  mixed  layer  are  separated  by  a  region  where  neither  U*  nor  z^ 
are  Important  and  free  convection  scaling  might  apply,  with  z,  g/Tv  and 
H0/p0cp  as  the  governing  variables.  However  over  the  oceans  the  surface 
heat  flux  is  much  smaller  than  over  land,  consequently  the  flow  in  the 
large  region  between  the  surface  layer  and  mixed  or  convective  layer  is 
governed  not  only  by  those  parameters  but  also  U*  and  z^  as  well.  Using 
the  U*  as  scaling  velocity,  the  dimensionless  turbulence  parameters 
should  then  depend  on  two  independent  parameters,  z/z^  and  z^/L.  As  one 
examines  the  literature,  turbulence  data  rom  the  A.B.L.  are  usually 
presented  in  both  surface  layer  and  mixed  layer  scaling  schemes. 

7.1  Mean  Mind  Profiles 

7.1.1  Observations  Above  the  Surface  Layer 

The  flow  in  the  larger  part  of  the  ABL  above  the  surface  layer 
depends  not  only  on  the  elevation,  z,  and  the  surface  fluxes  of  momentum 
and  virtual  heat,  but  also  on  the  ABL  height  z^,  Corriolis  effects, 
atmospheric  conditions  above  the  ABL  and  entrainment.  Several 
similarity  schemes  have  been  introduced  but  none  of  them  have  included 
all  these  parameters.  Figures  2,  3,  4  and  5  show  mean  wind  profiles 
taken  under  different  stability  conditions.  Obviously  the  observations 
at  higher  elevations  above  the  surface  layer  can  only  be  made  by  balloon 
soundings  or  instrumented  aircraft,  and  therefore  are  not  made 
simultaneously.  The  velocity  profiles  shown  in  Figures  3  and  4  are 
averages  of  several  soundings.  Under  thermally  unstable  conditions, 
trade-wind  velocity  profiles  show  velocity  maxima  to  occur  well  within 


the  ABL  (Figs.  2  and  3),  while  profiles  measured  at  larger  latitudes  do 
not  seem  to  exhibit  these  maxima  (Fig.  4).  In  the  presence  of  low-level 
Inversions  with  warmer  air  blowing  over  cooler  water,  mean  wind  profiles 
show  well  defined  maxima  at  very  low  elevations  (Figs.  5,  6).  Under 
these  conditions,  as  the  turbulence  is  suppressed,  Coriolis  effects 
become  Important  and  conditions  vary  rapidly  as  the  observations  of 
Figure  5  indicate. 

From  these  observations  one  can  conclude  that  a  mixed  layer  with 
near-uniform  velocity,  as  Is  often  observed  over  land  surface,  is  seldom 
present  over  the  ocean,  and  Instead  wind  shears  often  exist  throughout 
the  ABL.  It  has  been  suggested  that  the  effect  of  barocllnlclty 
(horizontal  temperature  gradient),  not  Included  as  a  governing  parameter 
In  the  mixed  layer  similarity  scheme,  may  be  responsible  for  the 
observed  mean  wind  maxima  In  the  ABL  (28). 

7.1.2  Observations  In  the  Surface  Layer 

In  the  lowest  part  of  the  ABL  (Surface  Layer)  the  flow  is 
characterized  by  a  rapidly  decreasing  mean  velocity  and  turbulence 
toward  the  earth's  surface.  Near  the  surface  the  mechanically  produced 
turbulence  dominates  the  flow  while  thermal  stratification  has  an 
important  effect  on  the  turbulence  only  in  extreme  situations. 

Generally  speaking,  the  surface  layer  flow  over  land  surfaces  is 
affected  by  the  large  variations  in  surface  heat  flux  during  day  and 
night.  This  heat  flux  (upward  during  the  day)  is  responsible  for  the 
generation  of  convective  turbulence,  but  creates  stable  conditions 
(downward  heat  flux)  during  night  when  turbulence  is  suppressed.  Over 
the  ocean  these  diurnal  variations  are  much  less  pronounced  and  the  heat 


flux  Is  generally  much  smaller  than  over  land.  Also  the  heat  flux  over 
the  ocean  is  dominated  by  the  flux  of  latent  heat  (water  vapor)  rather 
than  sensible  heat. 

In  addition,  it  has  been  observed  that  the  ocean  surface  is 
aerodynamical ly  speaking  extremely  smooth.  Even  under  strong  wind 
conditions  U  >  12  m/s  in  the  presence  of  large  waves,  the  roughness 
length  zQ  is  less  than  1mm  [29|.  For  flow  over  a  land  surface  this 
value  of  Zg  would  be  equivalent  to  that  of  flat  desert  terrain  [11]. 

The  influences  of  heat  flux  and  roughness  on  the  mean  velocity  profile 
are  shown  schematically  in  Figure  7. 

Close  to  the  water  surface  the  flow  must  be  influenced  by  the  waves 
themselves.  The  extent  of  the  influence  of  the  waves  on  the  air  flow 
directly  above  is  not  exactly  established  but  must  include  among  other 
things,  the  wave  amplitude,  which  itself  is  governed  by  the  wind  speed. 

Similarity  theory  applicable  to  the  surface  layer  requires  the 
assumptions  of  steady  state,  horizontal  homogeneity  and  the  fluxes  to  be 
inependent  with  height.  Under  the  condition  of  neutral  stability  in  the 
absence  of  any  vertical  heat  flux,  the  theory  leads  to  the  well  known 
logarithmic  wind  profile. 

U*  =  £  in(z/z0>  (24) 

where,  Zg,  the  roughness  length  must  be  characteristic  of  the  underlying 

roughness,  in  this  case  the  water  waves.  In  practice,  Zg  is  determined 

by  measuring  mean  velocity,  U,  at  several  elevations  and  extrapolating  a 

straight  line  through  a  plot  of  the  velocity  observations  versus  an  z 

down  to  the  axis  where  the  mean  velocity  is  zero.  However,  when  the 

extrapolation  is  performed  on  a  set  of  velocity  observations  obtained 


within  the  layer  where  the  waves  Influence  the  flow  directly,  the  values 
of  Zq  may  not  be  representative  of  the  actual  roughness. 

In  the  case  of  non-neutral  thermal  stratification,  the  logarithmic 
profile  equation  needs  to  be  modified  to 

h  =  k  <25) 

where  i|>  is  a  universal  function  based  on  the  stability  parameter  z/L 
(see  Section  6.5  of  Reference  11).  Many  of  the  values  of  zQ  presented 
in  the  literature  are  based  on  velocity  measurements  below  10m  and 
assuming  neutral  stability.  This  practice  leads  to  tremendous  scatter 
of  zQ.  For  instance,  values  of  zQ  obtained  from  more  than  1000  velocity 
profiles  acquired  over  the  Baltic  Sea  and  North  Sea  [29]  are  scattered 
between  1x10"^  mm  and  100  m.  Even  when  the  results  were  limited  to  only 
those  profiles  for  which  the  thermal  stratification  was  near  neutral, 
values  of  zQ  still  show  a  considerable  scatter  and  moreover  exhibit  a 
clear  dependence  on  either  the  shear  velocity  or  mean  velocity  at  a 
given  reference  height  [29,31]. 

In  a  later  paper  Krugermeyer  et  al.  [14]  show  that  close  to  the 
water  surface  (within  three  wave  heights),  wind  profiles  are  distorted 
to  lower  wind  speeds  in  comparison  to  velocity  profiles  over  a 
comparable  rigid  surface.  Consequently,  they  conclude  that  the 
determination  of  the  stress  and  roughnes  length  by  the  profile  method 
must  be  based  on  velocity  observations  above  this  layer  of  wave  height 
Influence.  Other  investigators,  including  Fujita  [32],  who  have  made 
multilevel  mean  velocity  measurements  close  to  the  sea  surface  observed 
appreciable  departures  from  the  logarithmic  profile  especially  under 
strong  wind  conditions  with  wave  heights  in  excess  of  2m.  Many  of  the 


profiles  show  a  single  velocity  maximum  below  10m,  while  under  strong 
wind  conditions  two  velocity  maxima  have  been  observed  (Fig.  8). 

Based  on  these  observations,  the  mean  wind  profile  in  the  ABL  over 
the  ocean  can  vary  in  extreme  fashion  depending  on  the  sea  state,  the 
virtual  heat  flux,  conditions  at  the  upper  edge  of  the  boundary  layer  as 
well  as  the  general  atmospheric  conditions  in  the  free  atmosphere. 
Consequently  no  simple  similarity  model  is  available  at  the  present  time 
to  describe  the  ABL  flow,  and  be  able  to  predict  the  mean  velocity 
profile  with  any  certainty. 

7.2  Velocity  Variances  and  Momentum  Flux 

In  general  the  atmospheric  motions  are  categorized  as  macroscale, 
mesoscale  and  microscale,  where  the  latter  is  considered  the  turbulent 
part  of  the  flow.  Obviously  it  is  the  turbulence  in  the  ABL  that 
affects  the  operation  of  low-velocity,  low-altitude  aircraft.  Since  the 
turbulence  is  a  random  phenomenon  it  is  best  studied  by  variances, 
spectra  and  scales.  In  the  ABL  the  mechanically  produced  (small  scale) 
turbulence  dominates  the  flow  near  the  surface  but  decreases  rapidly 
with  height  as  the  large-scale  convective  turbulence  becomes  more 
important  (Fig.  9). 

Above  the  ocean  the  convective  turbulence  created  through  buoyancy 
effects  is  relatively  weak.  Not  withstanding,  a  well-mixed  boundary 
layer  can  exist  as  a  result  of  large-scale  motions  such  as  roll  vortices 
which  are  primarily  responsible  for  the  upward  transport  of  turbulent 
and  moist  air.  As  this  moist  air  is  moved  upward  it  cools  and  will 
eventually  condense  and  form  the  typical  cloud  streets  associated  with 
these  roll  vortices  (Fig.  10).  These  large-scale  structures  including 
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gravity  waves,  which  are  frequently  observed  in  the  ABL  under  thermally 
stable  conditions  (34],  cannot  be  considered  turbulence.  These  low- 
frequency  motions  are  nearly  two-dimensional  affecting  the  vertical 
turbulence  component  very  little.  Other  large-scale  phenomena  include 
thunderstorms  and  squalls,  etc.,  which  are  usually  of  short  duration  and 
characterized  by  a  sudden  increase  in  windspeed  and  turbulence  followed 
by  a  sudden  decrease  as  shown  in  Figure  11. 

The  turbulence  characteristics  of  the  ABL  discussed  in  the 
following  sections  will  apply  primarily  for  stationary  f low-observations 
or  for  those  observations  where  the  large-scale  fluctuations  are  removed 
from  the  data  through  some  form  of  high-pass  filtering. 

7.2.1  Observations  Above  the  Surface  Layer 

In  Figures  12  and  13  the  standard  deviations  of  the  velocity 

component  and  the  turbulent  stress  are  presented  for  two  different 

stability  cases.  The  values  of  o  (a  =  u,v  and  w)  have  also  been 

a 

normalized  by  the  friction  velocity,  U*,  conforming  to  surface-layer 
similarity  and  by  the  convective  velocity,  W*,  for  comparison  in  the 
convective  outer  layer  of  the  ABL.  The  observations  in  the  ABL  below 
the  cloud  base  are  quite  similar  for  the  two  cases  but  considerable 
differences  for  both  the  standard  deviations  and  the  stress  exist  above 
the  ABL.  Near  the  surface,  ou  increases  rapidly  as  the  surface  is 
approached,  o and  ow  increase  also  but  much  less  rapidly.  For  the  near 
neutral  case,  the  surface  values  of  the  ratios  ou/U*,  av/U*  and  c^/U* 
are  2.3,  1.6  and  1.25  respectivey,  in  reasonable  accordance  with  the 
values  observed  over  flat  terrain  [11].  For  the  more  unstable  case 
(Fig.  13)  with  increased  cloudiness,  the  surface  values  of  o  /U*  are 
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3.8,  2.3  and  1.8  respectively  which  are  systematically  higher  than  those 
for  the  near  neutral  case.  Above  100m  the  standard  deviations  decrease 
but  become  nearly  constant  in  the  upper  1/3  of  the  ABL,  where  the  values 
of  a  /W*  are  0.95,  0.55  and  0.4  respectively  for  the  near  neutral  case. 

a 

In  both  the  near-neutral  and  the  more  unstable  case  the  turbulence 
components  do  not  vanish  at  the  edge  of  the  boundary  layer.  In  the 
near-neutral  case  the  standard  deviations  of  the  3  velocity  components 
just  below  the  cloud  base  are  approximately  70%  of  their  respective 
surface  values.  On  the  other  hand  in  the  more  unstable  case  with  the 
presence  of  clouds,  the  standard  deviations  of  the  three  velocity 
components  are  approximately  equal  at  the  edge  of  the  boundary  layer. 

Figures  12  and  13  also  show  the  distribution  of  the  turbulent 
stress  -uw  in  the  ABL  which  in  both  cases  varies  linearly  with  height 
decreasing  to  about  zero  at  the  edge  of  the  boundary  layer. 

Measurements  made  during  the  GARP  Atlantic  Tropical  Experiment  in 
the  summer  of  1974  reveal  that  the  distribution  of  ow  is  quite  dependent 
on  stability  (Fig.  14),  and  that  the  normalization  with  the  convective 
velocity  certainly  does  not  collapse  the  data.  The  shape  of  the  aw/W* 
distribution  resembles  that  of  Figure  12,  although  the  values  of  the 
bulk  stability  parameter,  z^/L,  except  for  day  218  are  all  considerably 
smaller  than  -10.  The  observations  at  the  lowest  level  (approximately  z 
=  20m)  show  average  values  of  z/L  of  -0.5,  -1.5,  -1.5  and  -1.2 
respectively  for  the  four  days,  with  average  values  of  °W/U*  for  these 
days  all  in  excess  of  2.0. 

Measurements  of  the  standard  deviations  and  stress  made  over  the 
North  Sea  on  February  26,  1975  [36]  show  remarkable  similarity  with  the 
trade  wind  observations  (Fig.  15).  Later  observations  over  the  North 
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Sea  around  Great  Britain  (25|  under  a  variety  of  meteorological 
conditions  are  shown  in  Figure  16.  This  multiple  data  set  was 
subdivided  in  two  categories,  class  A  with  z^/L  <  -1.9  (average  -3.9,  46 
runs)  and  the  near  neutral  class  for  which  z^/L  >  -1.5  (average  -0.9, 
twenty  eight  runs).  The  distributions  of  a^/U*  shown  in  Figure  16  for 
the  two  classes  show  clearly  the  variation  with  stability,  especially 
for  the  lateral  velocity  component.  The  standard  deviation  of  the 
vertical  component  ow/U*  seems  to  be  independent  with  stability  in  the 
surface  layer  (z/z^  <  0.1),  but  shows  considerable  variation  with 
stability  above  the  surface  layer.  The  maximum  value  of  o^/U*  at 
z/z^  *  0.5  and  the  distribution  of  aj U*  are  similar  to  the  results 
shown  in  Figure  14.  The  lateral  component,  ov/U*,  depends  rather 
strongly  on  the  stability  parameter  z^/L  for  all  heights.  Figure  16 
also  shows  for  comparison  purposes  the  average  values  of  aa/U*  measured 
in  the  surface  layer  over  water  by  eight  different  investigators.  The 
stress  distribution  seems  relatively  independent  with  stability  as  the 
results  of  Figure  16  indicate  and  also  decreases  uniformly  with  height 
as  was  the  case  with  previous  observations  (Figs.  12,  13,  15). 

7.2.2  Observations  in  the  Surface  Layer 

Observations  of  the  atmospheric  flow  in  the  surface  layer  over 
water  are  much  more  numerous  and  are  available  for  a  large  variety  of 
atmospheric  and  stability  conditions.  Since  the  detailed  knowledge  of 
the  mean  and  turbulent  flow  in  the  lower  50  m  is  required  for  the 
operation  of  low-velocity  aircraft  near  the  ocean  surface,  the  remainder 
of  the  analyzed  data  is  primarily  based  on  observations  made  below  this 


Table  1  provides  the  reader  with  a  summary  of  all  the  references 
used  In  this  report,  from  which  turbulence  and  mean  velocity  data  have 
been  obtained.  Including  specific  Information  as  to  the  location. 
Instrumentation  and  platform  used  during  the  experiment.  Also  the 
height  of  the  observations  above  the  water  surface  as  well  as  the  range 
of  stability  conditions  and  the  mean  velocity  range  and  the  number  of 
observations  analyzed  are  listed  In  this  table.  Average  values  of  the 
variance  ratios  a  /U  and/or  a  /U*  were  taken  directly  from  the  text  for 

a  a 

eight  different  references.  However,  for  the  majority  of  the  listed 
references,  the  number  of  observations  of  the  mean  velocity,  shear 
velocity  and  the  standard  deviation  of  the  turbulence  components  are 
presented  In  this  table. 

The  majority  of  the  low-level  observations  were  made  below  15  m,  of 
which  some  were  made  above  shallow  water  on  fixed  masts.  However,  most 
of  the  observations  over  the  open  ocean  and  below  15  m  were  either 
acquired  from  the  stable  platform  operated  by  the  Bedford  Institute  of 
Oceanography  or  from  the  resarch  platform  F.L.I.P  (Floating  Instrument 
Platform)  operated  by  the  Marine  Physical  Laboratory  of  the  Scripps 
Institute  of  Oceanography  (Fig.  17).  For  higher  elevations  the  results 
are  primarily  based  on  observations  made  with  instrumented  aircraft  and 
on  observations  made  with  F.L.I.P.  at  the  30  m  level  and  on  the 
observations  acquired  from  the  miorometeorological  tower  on  the  beach  at 
Wallops  Island,  Virginia. 

7.2.2. 1  General  Turbulence  Observations 

Because  of  the  variations  in  thermal  stability,  atmospheric 
conditions,  instrumentation  and  data  analysis,  discrepancies  between  the 


different  sets  of  data  are  unavoidable.  Average  values  of  o^/U*  and 
oa/U  are  presented  for  each  elevation  In  Tables  2  through  7  and  Figs.  18 
through  23.  Also  values  of  the  constants  Aq  and  m^  obtained  from  least 
squares  curve  fitting  of  oa  versus  U*  and  U  respectively,  are  presented 
In  these  tables.  Where  ever  possible,  values  of  the  friction  velocity, 
U*,  obtained  with  the  direct  measurement  of  uw  were  compared  with  values 
of  the  friction  velocity  as  obtained  with  either  the  mean  velocity 
profiles  or  with  the  inertial  dissipation  technique.  However  based  on 
the  discussion  in  Section  6.2  of  this  report,  and  since  most  of  the 
observations  of  the  friction  velocity  were  obtained  with  the  use  of  the 
direct  measurement  method  of  uw,  preference  was  given  in  this  report  to 
those  values  of  the  friction  velocity  obtained  with  this  method. 

For  several  data  sets  enough  observations  were  available  to  study 
the  effect  of  stability  and  wind  speed  on  the  measured  turbulence 
quantities.  For  the  analysis  of  the  turbulence  data  only  local  values 
of  U*  and  U  were  used  for  normalization,  no  attempt  was  made  to  correct 
these  velocities  to  the  often  used  reference  height  of  10  m.  The 
average  values  of  aj\l  and  aa/U*  as  listed  in  Tables  2  through  7  are 
also  plotted  in  Figs.  18  through  23.  The  plots  of  a  /U  show 
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considerable  scatter  primarily  due  to  variation  in  stability.  The  solid 
lines  in  Figs.  18,  19  and  20  which  were  drawn  through  the  data  should 
provide  the  reader  with  reasonably  accurate  estimates  of  the  average 
turbulence  intensities  under  near  neutral  conditions.  The  scatter  of 
the  results  is  also  enhanced  by  the  fact  that  above  a  certain  critical 
mean  velocity,  turbulence  intensities  increase  with  wind  speed  as  a 
result  of  Increased  wave  height  and  aerodynamic  roughness  under  strong 
wind  conditions.  However  the  data  presented  in  Figs.  18,  19  and  20  are 


averaged  turbulence  intensities  for  which  the  mean  velocities  fall 
generally  below  this  critical  wind  speed.  Consequently,  the  scatter  of 
the  results  in  these  figures  must  be  primarily  due  to  variation  in 
thermal  stability.  Below  a  height  of  20  m  above  the  water  surface  the 
average  turbulence  intensities,  (0,,/U).^,  vary  between  6  and  12  %, 

U  avy 

(CTv/U)avg  between  5  and  12  %  and  (°w/u)avg  between  3  and  7  %. 

The  distribution  of  the  average  values  of  o  /U*  are  shown  in  Figs. 

a 

21,  22  and  23.  Especially  the  scatter  for  the  horizontal  velocity 

components  is  appreciably,  while  on  the  other  hand  the  scatter  of 

(a.  /U*)  „  is  much  more  limited.  It  Is  of  note  to  observe  in  these 
w  avg 

figures  that  the  average  values  of  c^/U*  over  flat  and  uniform  terrain 
(°U/U*  =  2.39,  av/U*  =  1.29,  o^/U*  *  1.25,  Ref.  11)  are  systematically 
lower  than  those  observed  over  water. 


7.2 .2.2  Variation  of  o a  and  U*  with  Mean  Velocity,  U 

Figures  24  through  28  show  clearly  how  the  slope  of  the  plots 
of  (a  =  u,  v  and  w)  and  U*  versus  U  change  at  a  mean  velocity  of 
approximately  10  m/s.  For  mean  velocities  below  this  level  the  straight 
line  fitted  through  the  data  goes  approximately  through  the  origin. 
However,  for  velocities  above  10  m/s  lines  of  best  fit  through  the  data 
intersect  the  ordinate  (a^-axis)  well  below  the  origin.  Consequently, 
the  turbulence  Intensities,  o  /U,  and  the  ratio  U*/U  and  therefore  the 

a 

drag  coefficient,  U*^/U^  are  velocity  dependent  above  some  critical  mean 
velocity,  which  Is  approximately  10  m/s  for  observations  at  a  height  of 
10  m  above  the  water  surface. 


The  data  from  References  17,  47,  48,  49,  50,  51,  53  and  6?  all  show 


a  similar  behavior.  Between  wind  speeds  of  10  m/s  and  20  m/s  the 


average  Increase  for  ou/U  is  38%  while  increases  of  35%,  28%  and  33%  are 
observed  for  o^/U,  ow/U  and  U*/U  respectively.  The  data  from  References 
17  and  48  both  taken  on  Sable  Island  show  an  excessive  increase  in  oy/U 

of  88%  In  the  same  wind  speed  range  and  therefore  have  not  been  included 

In  the  average  value  of  the  other  data  sets  of  35%.  Similarly,  the  drag 

coefficient  CD  =  U*2/U2  Increases  by  an  average  of  75%  based  on  the  data 

of  References  17,  58  and  53.  SethuRaman  (49)  attributed  the  abrupt 
increase  of  the  velocity  variance  with  mean  wind  speed  above  the 
critical  wind  speed  to  the  presence  of  large  roll  vortices,  since  only 
the  variances  of  the  horizontal  velocity  components  were  affected. 
However,  careful  1  analysis  of  the  observations  from  other  sources  show 
that  <jw  versus  U-plots  also  exhibit  abrupt  changes  in  slope  at  a  mean 
velocity  of  approximately  10  m/s.  Moreover,  Identical  variations  can  be 
observed  from  the  data  acquired  over  Lake  Flevo  in  the  Netherlands  [471, 
a  location  where  roll  vortices  are  unlikely  to  exist.  In  addition, 
similar  abrupt  changes  in  the  covariance  uw,  and  therefore  the  friction 
velocity  as  defined  by  U*  =  V-uw,  exist  for  the  observations  from 
References  17,  48  (Fig.  27)  and  53.  These  observations  seem  to  furnish 
more  evidence  for  the  increased-roughness  hypthesis  than  for  the 
existence  of  roll  vortices  as  an  explanation  for  this  phenomenon.  In 
Table  8  a  summary  is  presented  of  the  variation  of  aQ  with  mean  velocity 
below  and  above  the  critical  wind  speed,  while  in  Table  9  the  average 
turbulence  intensities  in  the  high-velocity  range  in  excess  of  10  m/s 
are  listed.  The  results  clearly  Indicate  the  increase  in  intensity  with 
mean  velocity  for  the  horizontal  components  and  also  a  slight  increase 
in  intensity  for  the  vertical  component. 


Figures  28  through  31  show  the  variation  of  <jq  and  U*  with  mean 
velocity  for  three  different  stability  categories  based  on  the  results 
listed  In  Reference  53.  Since  most  of  the  results  are  for  strong  wind 
conditions  (U  >  10  m/s)  the  lines  of  best  fit  do  not  go  through  the 
origin  but  well  below  it.  In  addition,  the  plots  exhibit  a  difinite 
stability  dependence  which  Is  not  negligible  especially  in  case  of  the 
two  horizontal  components.  The  effect  of  stability  in  the  case  of 
the  ow  and  U*  versus  U  plots  [Figs.  30,  31]  is  much  less  pronounced  and 
consequently  the  scatter  of  the  data  is  much  less  than  for  the  lateral 
and  streamwlse  components.  Based  on  these  observations  It  must  be 
concluded  that  the  turbulence  Intensities,  a  /U,  as  well  as  the  ratio 
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U*/U  (and  therefore  also  the  drag  coefficient  U*^/U^)  vary  in  varying 
degrees  with  mean  wind  speed  (U  >  10  m/s)  as  well  as  thermal  stability. 

The  thermal  stability  dependence  of  the  turbulence  intensities  and 
the  velocity  ratio,  U*/U,  can  also  be  ascertained  from  the  results 
plotted  In  Figs.  32  through  35.  Here  the  data  of  Reference  53  were 
classified  according  to  nine  stability  groups,  z/L  <  -0.10,  -0.10  <  z/L 

<  -0.05,  -0.05  <  z/L  <  -0.025,  -0.025  <  z/L  <  -0.01,  | z/L |  <  -0.01,  0.01 

<  z/L  <  0.025,  0.025  <  z/L  <  0.05,  0.05  <  z/L  <  0.10  and  z/L  >  0.10. 

For  each  stability  category  the  average  stability,  (z/L) avg  and  the 
straight  line  of  best  fit  of  oq  and  U*  versus  U  were  determined.  From 
the  latter  the  turbulence  intensities  a  /U  and  the  ratio  U*/U  were 
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calculated  for  given  mean  wind  speeds  between  10  m/s  and  20  m/s.  These 
results  generally  exhibit  an  increase  with  mean  velocity  and  a  decrease 
with  stability.  It  is  surprising  to  note  that  aw/U  Is  relatively 
insensitive  to  stability  but  increases  with  mean  wind  speed.  Also  of 
note  is  the  rather  eratic  behavior  of  U*/U  with  stability.  The  results 


plotted  In  the  Figs.  32  through  35  show  the  following  variation  due  to 
velocity  and  stability,  7.5%  <  ou/U  <  12%,  6%  <  oy/U  <  11%,  4%  <  c^/U  < 
6%  and  0.032  <  U*/U  <  0.043. 

A  different  way  of  analyzing  the  data  is  to  group  the  turbulence 

Intensities  and  the  shear  ratio  U*/U  in  velocity  categories  of  1  m/s 

intervals.  As  an  example,  all  the  available  data  between  the  mean  wind 

speeds  of  11.0  and  12.0  m/s  are  plotted  versus  z/L  in  Fig.  36.  These 

results  show  a  clear  stability  dependence  for  o^/U  and  av/U,  both 

increasing  with  decreasing  stability,  while  the  ratios  ow/U  and  U*/U  are 

relatively  Insensitive  with  stability  over  the  stability  range  shown. 

Similar  plots  for  the  other  velocity  categories  between  11.0  and  20.0 

m/s  show  identical  behavior,  although  for  mean  wind  speeds  approaching 

20  m/s,  the  stability  range  for  the  available  data  is  much  narrower  and 

the  stability  dependence  Is  less  obvious.  For  mean  wind  speeds  above  20 

m/s  the  results  show  actually  a  slight  decrease  with  decreasing 

stability  over  the  stability  range  j z/L |  <  0.05  (Fig.  37).  The  general 

stability  dependence  can  be  established  by  the  fitting  of  straight  lines 

to  the  data  in  each  velocity  category.  The  values  of  a  /U  and  U*/U  for 

a 

zero  stability  (z/L  =  0)  are  plotted  against  the  average  velocity  for 

each  velocity  category  (Fig.  38).  Average  values  of  a  /U  for  four 

a 

different  mean  wind  speeds  and  three  stability  values,  z/L  =  0.2,  z/L  = 
0.0  and  z/L  =  +0.2,  based  on  the  data  of  Reference  53  are  also  presented 
in  Table  10. 

7.2.2. 3  Variation  of  o  with  U* 

a 

Special  emphasis  Is  placed  in  this  article  on  the  behavior  of 
velocity  ratios  a  /U*  and  on  the  coefficients  A  in  the  expression 


a  =  A  U*,  with  stability,  specifically  under  near  neutral  conditions, 
a  a 

Monin-Obukhov  similarity  theory  predicts  that  in  the  surface  layer 
velocity  ratios  oa/U*  are  functions  of  the  stability  parameter  z//L, 
independent  of  the  roughness  of  the  underlying  surface  and  wind  speed. 
Based  on  observations  over  uniform  terrain  and  primarily  under 
convective  conditions,  the  Monin-Obukhov  theory  is  adequate  to  explain 
the  behavior  of  the  vertical  velocity  component.  However  in  Reference 
[61]  it  is  pointed  out  that  ratios  of  the  standard  deviations  of 
horizontal  velocity  components  to  the  friction  velocity,  U*,  depend  on 
z^/L,  where  zn-  is  the  height  of  the  lowest  Inversion.  However,  for  most 
of  the  available  data  no  values  for  z^  are  given  and  approximately  40% 
of  the  data  have  no  values  for  z/L  listed.  The  ratios  o  /U*  or  the 

a 

coefficients  should  be  constants  for  near-neutral  thermal  stability, 
usually  assumed  to  occur  under  sufficiently  strong  wind  conditions  when 
z/L  and  z^/L  both  approach  zero  and  convection  in  the  surface  layer  is 
relatively  unimportant. 

First  we  will  investigate  if  the  velocity  ratios  a  /U*  vary  with 

a 

wind  speed  as  was  the  case  for  the  velocity  ratios  o  /U  and  U*/U.  The 

a 

plots  of  oa  versus  U*  did  not  indicate  any  change  from  a  linear 
relationship  and  always  extrapolated  towards  the  origin.  In  order  to 
study  the  effect  of  wind  speed  more  closely,  all  available  values  of 
a  /U*  were  divided  in  velocity  categories  of  1  m/s  interval,  and  a 

a 

straight  line  relation  aa/U*  =  m(z/L)  +  b  was  fitted  to  the  data  from 
which  values  of  oa/U*  at  neutral  stability,  z/L  =  0,  for  each  velocity 
category  were  obtained.  The  results  of  this  analysis  are  shown  in  Fig. 
39,  clearly  indicating  the  independency  of  the  ratios  a  /U*  with  wind 
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speed.  The  average  values  of  these  results  which  are  cu/U*  = 

2.71,  ov/U*  =  2.22  and  ow/U*  =  1.39  fit  the  results  over  the  entire 
velocity  range  quite  well.  Comparable  values  under  near-neutral 
conditions  obtained  from  observations  at  a  number  of  locations  over 
uniform  land  terrain  are  ou/U*  =  2.39,  ov/U*  =  1.92  and  ow/U*  =1.25 
[11].  Comparison  of  these  results  show  clearly  that  the  average  values 
obtained  above  water  surfaces  are  systematically  higher  by  approximately 
15*. 

The  stability  dependence  of  the  turbulence  ratios  oa/U*  re  shown  In 

Fig.  40,  where  the  data  of  Reference  53  were  classified  according  to  the 

same  9  stability  groups  as  before  in  article  7. 2. 2. 2.  The  average  value 

of  o^/U*  for  each  group  was  then  plotted  versus  the  average  stability. 

The  results  show  that  au/U*  and  oy/U*  both  Increase  under  increasing 

2 

thermally  unstable  conditions.  Both  o  ,/U*  and  U*  /oa,  are  quite 

w  u  w 

independent  with  stability.  Based  on  a  linear  regression  analysis  of 
these  results,  the  values  at  neutral  stability,  z/L  =  0,  are  au/U*  = 
2.71,  O  /U*  =  2.25,  o/U*  =  1.42  and  U*2/o,  o  =  0.266.  These  results 
compare  favorably  with  the  results  of  Fig.  39.  In  Fig.  41  all  the  data 
between  z  =  8  m  and  z  =  13.5  m  are  classified  in  discrete  stability 
categories  and  average  values  of  aa/U*  are  plotted  versus  average  values 
of  z/L.  For  this  set  of  data  the  stability  dependence  is  quite  similar 
to  that  observed  in  Fig.  40,  and  the  neutral  values  of  c^/U*  are  2.66, 
2.19  and  1.33  for  o  =  u,  v  and  w  respectively.  The  data  from  Reference 
51  which  were  acquired  over  the  open  ocean  and  which  are  summarized  in 
Fig.  42,  the  latter  taken  directly  from  this  reference,  show  identical 


behavior  and  neutral  values  of  a  /U*  of  2.63,  2.19  and  1.25  respectively 

CL 

for  a  =  u,  v  and  w. 

Figures  43,  44  and  45  show  plots  of  all  individual  observations  of 
a  / 11*  versus  stability  up  to  a  height  of  150  m.  Figure  43  shows  a 
strong  stability  dependence  for  the  data  below  a  height  of  20  m 
(normally  the  height  of  the  surface  layer)  and  in  the  stability  range 
| z/L |  <  0.5.  Similar  conclusions  can  be  drawn  for  the  results  in  Fig. 

44  where  the  0y/U*  observations  at  z  =  30  m  do  not  seem  to  fit  in  with 
the  observations  between  z  =  2.3  m  and  z  =  13.5  m.  Also  the  same 
observations  can  be  made  for  the  ow/U*  data  in  Fig.  45. 

The  plots  of  the  individual  observations  of  a  /U*  versus  z/L  and 

a 

restricted  to  the  surface  layer  (z  <  20  m)  In  Figs.  46,  47  and  48  show  a 
clear  stability  dependence  as  far  as  the  horizontal  turbulence 
components  is  concerned. 

The  plotted  data  as  well  as  the  linear  regression  line  clearly 
indicate  a  greater  stability  dependence  for  the  lateral  velocity 
component.  On  the  other  hand,  the  stability  dependence  of  the  vertical 
component  is  less  clear  and  actually  shows  an  increase  under  stable 
conditions  and  a  near  constant  value  under  unstable  conditions.  The 
results  published  by  DeLeonibus,  Fig.  13  in  Reference  63,  show  clearly  a 
quite  similar  behavior.  The  average  values  for  aa/U*  in  the  stability 
range  |z/L|  <  0.1  for  these  data  sets  are  2.64,  2.17  and  1.37 
respectively  for  a  =  u,  v  and  w.  In  Table  11  the  average  values  at 
neutral  conditions  and  obtained  by  different  methods  are  listed  together 
with  the  near-neutral  land  values  from  Panofsky  and  Dutton  [11 1 .  The 
values  of  o^/U*  for  data  acquired  over  water  are  12%,  15%  and  10%  higher 
as  compared  for  the  land  values  for  a  =  u,  v  and  w  respectively. 
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It  has  been  observed  that  in  the  surface  layer  over  land  surfaces 

the  turbulence  ratios  a^/ U*  and  oy/U*  are  increased  drastically  due  to 

the  presence  of  large-scale  fluctuations  of  the  order  in  the 

convective  boundary  layer  and/or  mesoscale  terrain  features  in  the  case 

of  complex  terrain,  and  that  in  both  cases  ow/U*  is  relatively 

unaffected  Ill].  Consequently  It  seems  reasonable  to  assume  that  the 

Increase  in  c»w/U*  over  water  is  primarily  due  to  the  travelling  waves 

which  certainly  affect  the  turbulent  flow  just  above  the  water 

surface.  Davidson  [13]  has  pointed  out  that  turbulence  ratios  a  /U*  not 

a 

only  vary  with  stability  but  also  with  the  wind-wave  coupling  parameter 
C/U*.  Here  C  is  the  phase  velocity  of  the  peak  energy-containing 
components  of  the  wave  spectrum.  Antonia  and  Chambers  [15]  observed 
that  Increases  in  o  /U*  are  the  combined  effect  of  increases  in  o  and  a 

a  a 

decrease  in  the  turbulence  stress,  provided  C/U*  >40.  While  for 
C/U*  <  40  values  of  ou/U*  and  crw/U*  are  In  reasonable  agreement  with 
those  measured  over  land  for  the  same  value  of  the  stability  parameter 
z/L.  During  the  Davidson  observations.  May  1969,  above  the  Atlantic 
Ocean  approximately  200  miles  northeast  of  Barbados,  conditions 
prevailed  for  which  C/U*  >  40.  However,  the  presented  values 
of  oa/U*  are  relatively  small  undoubtedly  as  a  result  of  the  short 
averaging  time  of  10.9  minutes.  Since  the  data  set  presented  in  this 
report  lacked  information  on  wave  speed,  one  can  only  assume  that  during 
these  experiments  conditions  prevailed  for  which  C/U*  >  40  and  observed 
values  of  a  /U  are  expected  to  be  in  excess  of  those  observed  over  land. 


7. 2. 2. 4  Turbulence  Integral  Scales 

Longitudinal  velocity  scales  are  usually  obtained  from  the  integral 
time  scale 

ao 

T  =  J  R(r)dT,  (26) 

0 

from  which  the  corresponding  length  scale  can  be  obtained  by  multiplying 
this  time  scale  by  the  mean  velocity  in  accordance  with  Taylor's  frozen 
wave  hypothesis.  This  kind  of  integral  scale  is  very  sensitive  to  the 
presence  of  low-frequency  velocity  fluctuations  and  consequently  on  the 
method  by  which  the  turbulence  is  separated  from  the  mean  velocity. 
Lateral  and  vertical  scales  can  be  obtained  by  integration  of  the  space- 
correlation  coefficients  with  lateral  and  vertical  separation  distances 
respectively.  In  order  to  obtain  reasonable  estimates  of  these  scales 
it  is  assumed  that  these  correlation  coefficients  decay  exponentially 
according  to 

R0(y')  =  expl-y'/Lj] 

and  (27) 

R  (z')  =  exp[-z'/Lz] 
a  a 

where  y'  and  z1  are  the  separation  distances  in  either  the  lateral  or 
vertical  direction.  Lateral  integral  scales  can  then  be  obtained  from  a 
least-squares  fit  of  expression  27  on  the  measured  correlation 
coefficients  for  various  separation  distances  y'  and  z'.  Some  of  these 
scales  have  been  measured  with  instrumented  towers  located  on  land  close 
to  the  water  [60,  64  and  65].  The  results  from  the  two  expermental 
sites  show  a  reasonable  agreement.  The  variation  in  the  longitudinal 


Integral  scales  measured  at  Wallops  Island  is  primarily  due  to  the 
presence  of  large  scale  low-frequency  velocity  fluctuations. 

An  alternate  method  for  obtaining  the  length  scale  of  the  energy- 
containing  eddies  is  by  determining  the  wavelength,  xm,  corresponding  to 
the  peak  of  the  logarithmic  spectrum  which  is  related  to  the  integral 
scale  as  follows  [66], 


L 


x 
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(28) 


where  nm  is  the  peak  cyclic  frequency  and  fm  the  peak  reduced 
frequency.  However  often,  primarily  under  convective  conditions  the 
logarithmic  spectral  functions  are  quite  flat  near  their  maximum  and  no 
distinct  peak  can  be  discerned  except  for  the  vertical  velocity 


spectra.  Typical  values  for  the  peak  reduced  frequency  for  the  w- 
spectra  are  fm  =  0.25  for  unstable  conditions  and  fm  =  0.5  for  stable 
conditions  [591.  With  expression  28  this  leads  to  values  of  L*  =  25  m 
and  L*  =  13  m  respectively  at  a  height  of  40  m  which  is  in  reasonable 
agreement  with  the  value  of  L*  =  32  m  from  Reference  65  (Table  12). 

So  far  our  discussion  of  integral  scales  has  been  based  on 
observations  made  from  stationary  platforms.  Spectra  measured  from 
moving  platforms  such  as  aircraft  flying  in  the  along-wind  direction  are 
related  to  those  obtained  from  stationary  platforms  with  the  aid  of 
Taylor's  frozen  wave  hypothesis.  The  along-wind  u  and  v  spectra  contain 
the  low-frequency  fluctuations  similar  to  the  spectra  obtained  from  a 
stationary  plat  form,  while  the  w  spectra  have  a  well  defined  peak  in 
both  cases  [36,67].  The  across-wind  u  and  v  spectra  lack  the  low- 


frequency  contribution  especially  under  near-neutral  conditions. 
Integral  scales  obtained  from  the  peak  wave  length  of  the  along-wind 


spectra  must  be  7'nterpreted  as  longitudinal  integral  scales 

like  L“  while  the  lateral  integral  scales  La  are  related  to  the  across- 
x  y 

wind  spectra.  The  Rsi/lts  presented  in  Table  12  show  that  the  along- 
wind  and  across-wind  integral  scales  from  aircraft  spectra  compare  quite 
well  with  the  integral  scales  obtained  from  the  correlation  coefficients 
from  an  array  of  stationary  wind  measuring  Instruments. 

8.  Summary  and  Conclusions 

8.1  General  Observations  and  Remarks 

The  purpose  of  this  report  has  been  to  provide  detailed  information 
of  the  atmospheric  flow  In  the  surface  layer  above  the  ocean,  say  below 
30  m.  This  information  is  required  as  input  for  the  design  and 
operation  of  remotely  piloted  aircraft  at  low  velocities  and  low 
altitudes  above  the  water  surface.  In  this  report  most  of  the  available 
experimental  data  are  reviewed  and  brought  together  in  order  to  provide 
this  information. 

Many  of  the  flow  observations  in  the  surface  layer  over  the  ocean 
suffer  from  experimental  error  and  lack  of  calibration  control  etc. 
Despite  continuous  improvements  of  the  measurement  techniques,  the  data 
handling  and  the  data  analysis,  one  must  expect  a  considerable  scatter 
of  the  results.  Some  of  this  scatter  is  experimental  in  nature,  some  is 
due  to  the  var 
Above  all  the 
time  and  space 

must  be  Interpreted  as  averages  based  on  a  limited  number  of 
observations.  What  the  probability  is  that  the  average  results  will 
actually  occur  at  any  given  time  and  at  any  instant  is  difficult  to 
ascertain.  In  order  to  get  some  idea  of  the  variability  of  the 


iable  conSI^ 


atmospheric  flow  must  be" 


t  the  sea-air  interface. 


.Heing  random  in 


and  consequently  the  results  presented  iTTTfi H -tgaart 
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observations,  standard  deviations  of  the  presented  data  have  also  been 
Included  In  the  tables. 

There  Is  a  number  of  important  differences  between  the  atmospheric 
flow  over  water  and  over  land.  The  most  obvious  difference  Is  that  of 
the  surface  which  Is  fixed  over  land  and  Is  In  motion  over  water, 
examplified  by  travelling  waves  generated  by  the  atmospheric  flow.  The 
Influence  of  the  atmospheric  flow  on  the  disturbed  water  surface  and 
vice  versa  the  influence  of  the  disturbed  water  surface  on  the  flow  is 
extremely  complex  and  not  well  understood.  Analysis  of  the  data  In  this 
report  show  that  mean  wind  profiles  and  turbulence  are  extremely 
variable  below  10  m.  Most  likely  as  a  results  of  wave  effects. 

8.2  Conclusions  for  the  Mean  Velocity  Profile 

The  variability  of  the  mean  flow  In  the  ABL  is  evident  of  the  data 
presented  in  Figures  2,  3a  and  4  and  5,  while  the  mean  flow  variability 
near  the  water  surface  is  well  demonstrated  in  Figures  6  and  8.  Despite 
a  stability  dependence,  the  variability  of  the  mean  velocity  profile  Is 
clearly  exemplified  by  the  scatter  of  the  profile  roughness,  zQ,  as 
presented  in  references  29  and  31.  In  both  cases  the  variability  of  zQ, 
based  on  many  hundreds  of  velocity  profiles  measured  near  the  water 
surface  below  z  *  10  m  is  several  orders  of  magnitude,  varying  between 
10"  ^  mi)  and  100  mm.  Krugermeyer  et  al.  [14]  showed  that  the  mean  wind 
profile  is  dlstorded  by  the  waves  to  heights  equivalent  to  three  wave 
heights.  Under  strong  wind  conditions  U  >  15  m/s,  the  significant  wave 
height  may  be  as  much  as  5  m. 
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8.3  Conclusions  for  the  Turbulence  Statistics 

Near  the  water  surface  the  turbulence  statistics  such  as  variances 
and  covariances  are  dependent  on  the  thermal  stability  of  the  atmosphere 
and  on  the  sea  state  which  itself  is  closely  related  to  the  prevailing 
wind  speed  and  the  origin  of  the  generated  wave  field.  The  significant 
conclusions  based  on  the  reference  survey  and  the  data  analyzed  in  this 
report  can  be  summarized  as  follows. 

8.3.1  Turbulence  Ratios  a  /U* 

a 

1.  Turbulence  ratios  oq/U*  in  the  surface  layer  over  water  (z  <  20  m) 
deviate  from  the  average  observations  over  land,  provided  the  phase 
velocity  C  of  the  energy  containing  waves  is  at  least  40  times  the 
friction  velocity  U*  [13]. 

2.  Under  the  above  conditions  increases  in  the  turbulence 

ratios  o  /U*  are  due  to  an  increase  in  o  and  a  decrease  in  U* 
a  a 

[151. 

3.  For  the  data  sets  used  in  this  report,  values  of  C/U*  are  generally 
not  available.  However,  based  on  the  Davidson  results  [13],  it 
seems  appropriate  to  assume  that  for  fully  developed  waves,  which 
generally  occur  over  the  open  ocean,  the  ratio  C/U*  is  generally 
larger  than  40.  Therefore  larger  values  of  °a/U*  can  be  expected 
in  the  surface  layer  over  water  than  over  land.  (Figs.  21,  22  and 
23). 

4.  The  turbulence  ratios  o  /U*  are  generally  independent  with  wind 

a 

speed  (Fig.  39). 

5.  The  stability  dependence  of  the  turbulence  ratios  oq/U*  is 
different  in  the  surface  layer  (z  <  20  m)  than  at  higher  elevations 
(z  >  20  m)  (Figs.  43,  44  and  45). 
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6.  Both  ou/U*  and  ov/U*  increase  appreciably  with  decreasing  stability 
(Figs.  40,  41  and  42). 

7.  In  the  surface  layer  (z  <  20  m)  the  linear  regressions  au/U*  =  2.64 
-  1.02  z/L  and  ay/U*  =  2.22  -  3.09  z/L  are  the  best  fit  to  the 
available  data  (Figs.  46  and  47). 

8.  In  the  surface  layer  (z  <  20  m)  the  turbulence  ratio  aw/U*  is 
relatively  Independent  with  stability  (Figs.  40,  41,  42  and  48). 

9.  In  the  surface  layer  (z  <  20  m)  and  for  near-neutral  stability 
conditions  average  values  for  cr^/U*  are  2.67,  2.20  and  1.38  for  o  = 
u,  v  and  w  respectively.  Average  over-land  values  are  2.39,  1.92 
and  1.25  (Table  11). 


8.3.2  Velocity  Ratios  aj U  and  U*/U 

1.  The  velocity  ratios  o  /U  and  U*/U  vary  with  height  (Figs.  18,  19 

a 

and  20),  with  wind  speed  (Figs.  24  through  31  and  38),  and  with 
thermal  stability  (Figs.  32  through  37). 

2.  The  standard  deviations  of  the  turbulence,  a  ,  and  the  friction 

a 

velocity,  U*  =  V-uw  vary  with  mean  wind  speed  in  different  ways 
depending  whether  the  mean  wind  speed  is  below  or  above  some 
critical  value  (Table  8). 

3.  This  change  in  mean  wind-speed  dependence  is  probably  the  result  of 
increased  roughness  which  occurs  under  strong  wind  conditions. 

4.  The  critical  value  for  U  varies  with  height  but  is  approximately  10 
m/s  at  z  =  10  m. 

5.  For  wind  speeds  below  the  critical  value,  oa  and  U*  vary  linearly 
with  U.  The  line  of  best  fit  generally  passes  the  origin  closely 
(Table  8). 


6.  For  wind  speeds  In  excess  of  the  critical  wind  speed,  i  and  u* 

3 

stm  vary  linearly  with  U.  However,  the  slope  of  the  linear 
regression  Is  Increased  and  deflr'tely  does  not  pass  through  the 
origin  (Figs.  24  through  31). 

7.  For  mean  velocities  below  the  critical  wind  speed  and  for  near 

neutral  thermal  stability  conditions,  values  for  the  turbulence 

Intensities,  a  /U,  can  be  otalned  from  Figures  18,  19  and  20. 
a 

Average  values  under  these  conditions  at  z  *  10  m  are  o^/u  * 

8X,  c>v/U  *  6 %  and  a  / U  -  4.2*  while  UVU  -  0.034. 

8.  For  near-neutral  conditions  and  for  wind  speeds  In  excess  of  the 
critical  wind  speeds,  values  for  the  turbulence 

Intensities  a  /U  vary  with  height  and  wind  speeds  (Table  9). 

a 

9.  Averaged  turbulence  intensities  at  z  -  13.5  m  for  four  specific 
values  of  mean  wind  speeds  (U  -  10,  1 3 ,  16  and  19  m/s)  In  excess  of 
the  critical  wind  speed  and  for  three  stability  conditions  U/i  * 
-0.2,  0.0,  +0.2)  are  prsented  In  Table  10. 

10.  At  z  -  13.5  m,  depending  on  stability  and  mean  wind 

speed,  a^/U  varies  between  6.6  and  '2.8*.  ?y/U  between  4.2  and 
10.9*,  and  o^/U  between  4.2  and  5.8*. 

8.3.3  Integral  Scales  and  Spectra 

1.  Integral  scales  obtained  from  auto  or  space  correlation 
coefficients  are  very  sensitive  to  t*e  presence  of  low-frequent y  or 
large-scale  velocity  fluctuations,  and  consequently  vary  great  , 
magnitude  (Table  12). 

2.  Integral  scales  obtained  from  the  -ave'ength  corresponding  to  the 
peak  of  the  logarithmic  velocity  spectra  compare  rasonably  wei' 
with  those  obtained  from  correlation  oeff'c'ents. 


3.  Over  water  the  measured  integral  scales  are  generally  several  times 
larger  than  those  observed  over  uniform  land  terrain. 

4.  The  shape  of  the  velocity  spectra  depends  primarily  on  stability 
and  In  the  case  of  aircraft  data  on  the  flight  direction  (along 
wind  or  across  wind)  (67]. 

5.  Independent  of  stability,  more  low-frequency  energy  (f  <  0.1)  Is 
observed  in  the  along-wlnd  measured  u  and  w  spectra  while  In  the 
high-frequency  range  f  >  0.1  more  energy  is  present  In  the  across- 
wind  measured  spectra  [67]. 

6.  The  shape  of  spectra  measured  in  the  along-wind  direction  agrees 
well  with  the  shape  of  velocity  spectra  measured  over  land  [67]. 

8.4  Final  Remarks 

The  meteorological  data  which  were  analyzed  in  this  report  were 
primarily  acquired  for  vertical  flux  studies  of  momentum,  heat  and  water 
vapor,  and  were  not  acquired  for  the  purpose  of  providing  design  input 
and  operation  criteria  for  low-flying  aircraft.  Despite  this  and  the 
limited  amount  of  data  available,  a  reasonable  description  of  the  mean 
and  turbulent  flow  In  the  surface  layer  over  the  ocean  has  been 
presented  in  this  report.  Notwithstanding  the  fact  that  the  presented 
results  do  not  provide  a  complete  picture  of  the  flow,  it  is  believed 
they  are  adequate  as  an  input  to  the  design  and  operation  of  remotely 
piloted  aircraft  flying  at  low-Reynolds  numbers  just  above  the  ocean. 
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FIGURE  2  Mean  Velocity  profiles  in  the  trade-wind  boundary  layer 

(a)  near  neutral  conditions,  (b)  intermediate  convection, 

(c)  stronger  convection.  U-west  component,  V-south  component 
Surface  wind  predominantly  from  east  (Reference  23). 
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FIGURE  3a  Averaged  mean  velocity  profiles  taken  at  the  small 
island  of  Anegada  north-east  of  Puerto  Rico 
(Reference  24). 


bottom  scales,  U(o)  m/s 


FIGURE  3b  Averaged  mean  velocity  profiles  taken  at  the  small 
island  of  Anegada  northeast  of  Puerto  Rico 
(Reference  24). 


irofiles  over  the  sea  around  Great 
-3.9  (more  unstable  category) 
is s  unstable)  (Reference  25) 


FIGURE  10  Schematic  diagram  showing  the  circulation  associated  with 
longitudinal  roll  vortices,  and  their  association  with 
cloud  streets.  Buoyant  eddies  (stippled)  drawn  into  the 
circulation  are  helped  to  the  condensation  level,  where 
they  can  form  clouds  (Reference  33). 
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FIGURE  12  Standard  deviations  a  (a  =  u,  v  and  w)  and  -uw  versus 
elevation  in  the  tradi  wind  boundary  layer  northeast  of 
Puerto  Rico.  (U*  =  0.45  m/s,  W*  =  0.73  m/s,  z/L  = 

-0.07  and  z^/L  =  -1.4  from  Reference  28). 
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FIGURE  13  Standard  deviations  o  (a  =  u,  v  and  w)  and  -uw  versus  elevation 
the  trade  wind  boundary  layer  northeast  of  Puerto  Rico  (U*  =0.3 
m/s,  z/L  =  -0.2  and  Z../L  =  -3.0  from  Reference  28). 
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FIGURE  14  Distribution  of  a  /W*  in  the  ABL  under  rather  unstable  conditions 
z^/L  <  -10  (Reference  35). 


curves  numerical  prediction. 


Distributions  of  (°a/u*)aVg  and  (-uw/U*( ;avg  in  the  ABL  for 

*  Class  A  •  Class  B  (see  text).  Horizontal  line  through  the 
data  points  indicate  a  measure  of  the  standard  deviation  of 
the  average  values.  The  crosses  at  z/zj  =  0  denote  average 
values  observed  over  water  (Reference  25). 


FIGURE  17  Diagram  of  FLIP  (Floating  Instrument  Platform)  showing  the 
experimental  arrangement.  Integrated  accelerometer  spectra 
and  the  recorded  ship's  heading  were  used  to  correct  the 
measured  velocity  components. 
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FIGURE  22  Average  value  of  c  /U*  versus  height  for  three  stability 
categories.  Average  land  value  1.92  [11]. 
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FIGURE  23  Average  values  of  a  /U*  versus  height  for  three  stability  categori 
The  large  value  of  a  /U*  at  Z  =  43  m  represents  only  a  single  data 
point.  Average  land  value  1.25  [111. 


ices 


Standard  deviation  a  versus  mean  velocity,  for  data  with  | z/L |  <  0.15  from 
References  17  and  48,  z  =  10  m. 
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FIGURE  26  Standard  deviation  a  versus  mean  velocity,  for  data  with  |z/L 
from  References  17  and  48,  z  =  10  m. 
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velocity  U*  versus  mean  velocity  for  all  observations  of  Reference 
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turbulence  intensity,  ov/U,  %,  versus  (z/L)  for  different  mean  veloci 
-velocity  range  and  obtained  from  the  linesaor  best  fit  through  the  data 
ility  category  (ov  =  m  U  +  b).  Data  from  Reference  53,  z  =  13.5  m. 


Sw  /  U 


a  a  a  • 

0  0  0  3 

_ _ _ — .  cn 

3  3  3  3 

SSV 
3  3  3  3 


*--*  03 

—1  ^ 
\  — J 

Ni  \ 

0) 

- -  fsl 

u  m  ^ 
c  CO  co 

<U 

•  r— 

s- 

CNJ  O 

V 

• 

<*- 

I  o 

Ol 

QC 

to  + 

CVJ 

c- 

•  to 

o  M 

C\J 

05 

II 

4-> 

o 

03 

o> 

-o 

•  >11 

05 

0) 

^  CD 

JC 

*  > 

4-> 

=0  <0 

s- 

o 

D  D 

4- 

o 

>»  I  -- 

4->  I  CM 


•r— 

i  * 

r- 

=> 

•i— 

JO 

CD 

05 

> 

-M 

03 

</5 

— 1  . 

-d 

4-> 

^  ' 

•r— 

5 

1 

</) 

<D 

| 

+-> 

•r- 

1 

+-> 

cn 

c 

f—  > 

03 

03 

3 

•  "S 

ctcm  _j 

N 

Q) 

II  N 

U 

C 

con 

<u 

f— 

03  • 

or 

^  o 

JO 

* 

i- 

ZD  + 

3 

4-> 

O  ^ 

<U 

-. _ *  • 

o> 

r— 

03 

s_ 

II 

<D 

> 

cn 

03 

> 

4-  • 

O  £  * 

=) 

cm  m 
O  •  S 
•I-  co  D 


CC 

=) 

CD 


(Su/Ux)  avg. 
(Sv/Ux)  avg. 
(Sw/Ux)  avg. 


«  / 


v 

« 

<  /  4 
/  << 
A  * 


1  / 
/ 

/  « 


'A  0 


t  &/Uft 

♦  fc/ux 

“  x  o^ux 

0  0^/71 

1 

1  I 

-  1  _ 

1 

-0.3 

1  1 
-0.1 

i 

0.1 

Z/L 

FIGURE  42  Average  values  of  a  / U*  versus  z/L  from  Reference  51,  Fig.  18 
Average  values  at  zero  stability  are  2.63,  2.19  and  1.25  for 
a  =  u,  v  and  w  respecti vely. 
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Variation  of  o  /U*  versus  z/L  for  all  observations 


Variation  of  a/U*  versus  z/L  for  all  observations 
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FIGURE  47  Variation  of  o  /U*  versus  z/L  in  the  surface  layer  (z  <  13.5  m).  Linear  regression 
line  a  /U*  =  2^.22  -  3.09  z/L. 


Variation  of  q ,/U*  versus  z/L  in  the  surface 


z,n  NuMber  of  Observations  Analyzed  Method  for  Obtaining 


z,*  Number  of  Observations  Analyzed  Method  for  Obtaining 


z,n  NuMber  of  Observations  Analyzed  Method  for  Obtaining 
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Stability:  U  Unstable 

N  Near  Neutral 
S  Stable 

n  Number  of  Observations 


TABLE  8 


Variation  of  ofl  with  Velocity  Below  and  Above  the  Critical  Wind  Speed 


Reference 

z,m 

Velocity 

Range 

m/s 

o  * 
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mu 

m  U  +  b 
a  a 

bu 

a  *  m'U 
a  a 

Body  of  Water 

47 

8.4 

10.1-14.5 

0.140 

-0.559 

- 

Lake  Flavo, 

The  Netherlands 

17,48 

10.0 

5.6-10.0 

10.0-21.2 

0.092 

0.145 

-0.098 

-0.584 

0.080 

Atlantic  Ocean, 

Sable  Island 

49 

10.0 

<12.0 

>12.0 

0.20 

-1.5 

0.08 

Atlantic  Ocean 

5  km  off  Long  Island 

50 

11.7 

6.7-11.5 

11.5-15.9 

0.091 

0.185 

-0.055 

-1.15 

0.086 

Atlantic  Ocean 

53 

13.5 

10.1-20.0 

0.134 

-0.545 

- 

Atlantic  Ocean 

17,48 

10 

5. 6-9. 8 
10.0-21.2 

Illy 

0.059 

0.173 

bv 

0.007 

-1.20 

0.060 

Atlantic  Ocean 

Sable  Island 

49 

10 

<12.0 

>12.0 

0.16 

-1.06 

0.08 

Atlantic  Ocean 

5  km  off  Long  Island 

50 

11.7 

6.7-11.5 

11.5-15.9 

0.076 

0.129 

-0.045 

-0.663 

0.072 

Atlantic  Ocean 

53 

13.5 

10.1-20.1 

0.138 

-0.826 

- 

Atlantic  Ocean 

17,48 

10 

5. 6-9.8 
10.0-21.2 

mw 

0.045 

0.068 

bw 

-0.024 

-0.18 

m' 

w 

0.042 

Atlantic  Ocean 

Sable  Island 

49 

10 

<12.0 

>12.0 

— 

0.046 

0.053 

Atlantic  Ocean 

5  km  off  Long  Island 

50 

11.7 

6.7-11.5 

11.5-15.0 

0.041 

0.041 

0.047 

0.120 

0.046 

0.050 

Atlantic  Ocean 

53 

13.5 

10.0-20.1 

0.072 

-0.299 

0.054 

Atlantic  Ocean 

TABLE  9 


Turbulence 

Intensities 
o  /U.  % 


Intensities,  o  /U,  for  Mean  Wind  Speeds  In  Excess  of 
the  Critical  Wind  Speed 

Mean  Velocity,  U,  m/s 

10  12  14  16  18  20  References 


TABLE  10 


Summary  of  the  turbulence  Intensities  In  the  high-velocity  range 
Results  are  based  on  curve  fits  of  the  data  presented  In 
Figures  33  through  35.  (z  *  13.5  m) 


Turbulence 

Component 


-0.2 


Stability,  Z/L 

0.0 


TABLE  11 


Summary  of  the  Turbulence  Ratios 

o  /U*  a  3  u,  v  and  w 
a 

- 

Source 

au/U* 

%/U* 

V  U* 

*  n 

Fig.  39 

2.71 

2.22 

1.39 

V, 

Fig.  40 

2.71 

2.25 

1.42 

Fig.  41 

2.66 

2.19 

1.33 

Ref.  51 
z  =  13  m 

2.63 

2.19 

1.25 

|z/L|  <  0.1 
z  <  20  m 

2.64 

2.17 

1.37 

Suggested 
averages 
over  water 

2.67 

2.20 

1.38 

Suggested 
averages 
over  1  and  [11] 

2.39 

1.92 

1.25 

Water:  ^  =  1.93 

—  *  1.59 

**"* 

°w 

°w 

V; 

°u 

Land:  —  =  1.91 

°w 

CTV 

—  =  1.54 
aw 

p* 

mm 


TABLE  12 

Turbulence  Integral  Scales  Measured  Over  Water 


Height  z,m 

Scale 

Magnitude,  m 

15 

*5 

60  -  150 

30 

$ 

90  -  240 

40 

204 

40 

$ 

32 

40 

13 

40 

50 

15 

*5 

20  -  40 

30 

$ 

20  -  50 

15 

$ 

10  -  20 

30 

$ 

20  -  30 

40* 

$ 

14  -  20 

40* 

L“ 

32 

40* 

35 

Integral  scales  obtained  from  spectra  measured  by  aircraft. 


Reference 

60 

60 

64 

65 
65 
64 
60 
60 
60 
60 
67 
67 
67 


END 
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